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Abstract: Currently there exists no known way to construct the Stress- Energy Tensor (T M „) of 
the produced medium in heavy ion collisions at strong coupling from purely theoretical grounds. 
In this paper, some steps are taken in that direction. In particular, the evolution of T^ v at strong 
coupling and at high energies is being studied for early proper times (r). This is achieved in the 
context of the AdS/CFT duality by constructing the evolution of the dual geometry in an AdSs 
background. Improving the earlier works in the literature, the two incident nuclei have an impact 
parameter b and a non-trivial transverse profile. The nuclear matter is modeled by two shock 
waves corresponding to a non-zero five dimensional bulk Stress- Energy Tensor Jmn- An analytic 
formula for T^ v at small r is derived and is used in order to calculate the momentum anisotropy and 
spatial eccentricity of the medium produced in the collision as a function of the ratio The result 
for eccentricity at intermediate t agrees qualitatively with the results obtained in the context of 
perturbation theory and by using hydrodynamic simulations. Finally, the problem of the negative 
energy density and its natural connection to the eikonal approximation is discussed. 
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1. Introduction 



There has been strong evidence from RHIC that the Quark-Gluon Plasma (QGP) created in heavy 
ion collisions goes through a strongly coupled phase [1] - [19]. In particular, hydrodynamic sim- 
ulations which describe the data successfully, require small shear viscosity compared to entropy 
density. This fact sends the message that very early after the collision the produced medium ex- 
hibits a strongly coupled behavior. Hydrodynamics simulations also use a short thermalization time 
of the order of 1 fm/c [15]. 

There have been theoretical efforts trying to obtain thermalization and more importantly 
isotropization at early times within the context of perturbative QCD ( pQCD) [20] - [26]. While 
research in this direction is still under way, there has not been a satisfactory explanation of the 
(early) thermalization in the framework of perturbation theory. Nonetheless, perturbation theory 
is not only a powerful tool for describing collisions but also an essential one. The reason for this 
is due to the fact that the produced system starts out being weakly coupled and described in the 
Color Glass Condensate (CGC) framework [27] - [47] and subsequently becomes strongly coupled 
and thermalizes; the process takes place over times of the order of 1 fm/c. As a result, we should 
look for alternative methods to pQCD in order to describe the intermediate stages of heavy ion 
collisions. 

Fortunately such a tool that may be applied to strongly coupled theories exists: The anti-de- 
Sitter / Conformal Field Theory correspondence ( AdS / CFT) was discovered in 1997 by Maldacena 
[48] and was soon quantified by Witten [49]. Since then, its validity has passed many physics 
tests while it has been widely applied in problems related to QCD (and recently to condensed 
matter physics). Many observables or interesting quantities in QCD have been calculated and many 
processes have been studied within the framework of AdS/CFT. The list includes applications of 
AdS/CFT to the heavy quark potential [50] - [57], to the evaluation of scattering amplitudes, in 
Deep Inelastic Scattering and to small-x physics [58] - [80] , to parton motion in the plasma, energy 
loss and radiation [81] - [87], heavy ion collisions, QGP and the hydrodynamic behavior of the 
produced medium [10] - [15] and [88] - [113] and many others. Regarding heavy ion collisions, the 
question of whether a thermal medium is produced may be attacked by looking at trapped surfaces 
in the gravity side. This has been done in [114] - [128] and references therein, and the method has 
been adapted from earlier work on general relativity. We stress that AdS/CFT does not deal with 
QCD but with M = 4 SYM theory with an SU(N) gauge group snd with N being large. In [129] 
one may find the values of several thermodynamic quantities computed by applying lattice methods 
with finite N at large coupling and comparisons with results obtained from AdS/CFT [130-133]. 
There are more differences between the two theories and certainly many similarities, especially 
at the deconfining limit. Reference [134] exposes these matters and reviews the progress in the 
field, starting from 't Hooft's large N field theories [135] and building up the discovery of the 
gauge/gravity correspondence and its applications to problems related to gauge theories and QCD. 

In this paper, we are interested in studying the problem of heavy ion collisions in the context of 
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the gauge-gravity correspondence. In particular, we are interested in calculating the Stress-Energy 
(SE) tensor of the produced medium and related quantities along with the way they evolve with 
time. This is achieved by finding the time evolution of the dual geometry in an AdS background 
since according to the duality, the metric tensor on the gravity side is dual to the (expectation 
value of the) SE tensor in the gauge theory side. We assume that the nuclei that collide have 
non-trivial transverse profiles and that the collision is not central but there is a nonzero impact 
parameter (b) involved 1 . To the best of our knowledge, this is the first time in the literature where 
in the framework of the AdS / CFT correspondence modeling of heavy ion collisions with transverse 
dynamics is taken into account in order to calculate the evolution of the SE tensor of the produced 
medium. The collision is assumed to take place at large 't Hooft coupling 2 and our analysis applies 
at early proper times in a way that is fully explained and quantified. 

The problem from the gravitational perspective is formulated as an initial value problem where 
the dual geometry is known in some time interval (in contrast to [15], [95,96]). This initial geometry 
is chosen as two shock waves moving towards each other and represents the nuclear matter to be 
collided. The evolution of the geometry, as determined by Einstein's equations, maps onto the 
evolution of the produced matter of the gauge theory in a way which is made precise by the 
AdS/CFT correspondence. The paper [136] quantifies the relation between the SE tensor of the 
gauge theory and the metric tensor on the gravity side. A discussion which elucidates this issue 
through particular examples may be found in [15] and [95]. Our work has partial overlap with [92]- 
[100] and may be regarded as an improvement of [93] which in turn was motivated by [15]. All these 
papers employ the results of [136] . We also note that the literature about shock waves in gravity is 
very rich; references [137] - [145] include only a subset of it. 



We organize our paper as follows. 



In section |2] we set up the problem, explain the precise way that the problem of heavy ion 
collisions is mapped onto a gravity problem in the context of AdS/CFT and state how the SE 
tensor of the gauge theory can be computed from the metric of the dual theory. 

In section ^ we talk about the transverse part of the nuclear profiles and how these may be 
modeled in terms of gravitational shock waves which correspond to a non-zero bulk (gravity) SE 
tensor. We argue that in order to have a well defined and non-trivial metric, meaning that it 
is finite both in the UV and the infrared (IR) in AdS bulk, a bulk source is necessary (compare 
with [146]). This is analogous to classical electrodynamics where the Laplace equation has nontrivial 
solutions which are finite at both the origin and at infinity only in the presence of a charge density. 



1 In fact, the impact parameter is required, for otherwise one has to face violation of conservation of the bulk 



stress-energy tensor (see (4.20)) and ultraviolet (UV) infinities (see subsection 7.4). 

2 Although this is not the case in real QCD as at the very early times, as has been already mentioned, pQCD 
seems to describe collisions satisfactorily. 
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In particular we derive the Green's function ( |3.9| ) of the scalar operator for AdSs which has the 
desired boundary conditions. 

In section |], we take into account back-reactions in the gravity description in order to have a 
conserved and therefore a well defined bulk SE tensor (see ( |4.15| ) and ( 4.19 )). Here, we take into 



account that the first bulk source moves in the gravitational field of the second and vice versa. This 
alters the trajectories of the bulk sources and induces a time evolution in the matter density in 
the bulk as figure [| suggests. In a diagrammatic approach, these corrections to the initial bulk SE 
tensor would correspond to the diagrams of figure ||. 

Section ^| includes Einstein's equations and the way these equations may be decoupled in order 
to be solved perturbatively. The program is to expand the field equations in the AdSs background 
in the presence of the two shock waves and solve them order by order in the strength of the shock 
waves. At each step (order of the expansion) we make sure that the bulk SE tensor is conserved 
(see section |j). The first correction to the metric corresponds to the diagrams of figures § and [5[ 
Obviously, according to the dictionary of the AdS/CFT, this correction to the metric maps onto 
a correction of the initial SE tensor of the colliding nuclear matter. Finding this correction is the 
ultimate goal of our work. 

In section ^| we exhibit the main details of the calculation and state the main result of this 
paper, that is the formula giving the gauge theory SE tensor of the matter produced in a collision. 
This is given by ( 6.18 ) while the details of the calculation are left for the two Appendices at the 



end of the paper. 

Finally, in section [7| we discuss the area of validity of our approximations and quantify the 
meaning of the fact that our approximation applies for early times. This is done in subsection |7.1| . 

In subsection |7.2| we estimate the energy density and the transverse pressures at certain kine- 
matical regions (see ( |7.13| ) and ( |7.19| )) where our analytic formula (see ( |6.18| ) and ( |7.25| )) simplifies 



to give compact expressions for the SE tensor 3 . In particular, for the regions J, J and JJJ, III of 
figure H we find that the energy density behaves as in ( [7.1 3|) and (|7.19|) respectively. The result for 



regions JJJ, III' exhibits a similar logarithmic dependence on proper time as the results of [147] 
- [149] although our result (at central rapidities) has an additional overall factor of proper time 
squared. In [147] - [149] the results are obtained in the context of pQCD. 

In subsection [773] we use the approximate results of previous subsection in order to derive 
approximate formulas for the momentum anisotropy and the spatial eccentricity of the produced 
medium. In the two particular limiting cases that we have investigated in subsections |7.2| , the 
transverse pressure components of the SE tensor turn out to be equal. Therefore we conclude that 
the momentum anisotropy (see ( |7.21| )) should be close to zero. In addition, for one of the two cases 



we estimate the spatial eccentricity as a function of the ratio ^ for intermediate proper times r (see 
( [7.240 and figure 0). Our result is qualitatively consistent with the results of [150,151] obtained 



3 Certainly our analytic formula is (more) accurate but unfortunately rather complicated: we find it useful to 
investigate particular limits of it. 
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from pQCD and the results of [1] and [152, 153] obtained by using hydrodynamic simulations. Our 
analytic formula for the SE tensor, equation ( p.!8| ), allows us to extend these results to a broader 
kinematical region than the two aforementioned limiting cases (of subsection [7.2]) ; this is left for a 
future project. 

In subsection [7.4| we summarize our conclusions: The evolution of the gauge theory SE tensor is 
constrained by causality 4 (see figure |^ and equation ( |7.25| )) in an intuitive way. Causality basically 
separates the evolution into three 5 areas which are shown in figure | (see also equation (|7.25|) ). 
Mathematically this is due to the presence of retarded propagators which induce ^-functions, and 
which in turn separate space-time according to figure |6|: the SE tensor at any given point of space 
evolves according to whether the propagation of the signal from the center of each nucleus (or 
both nuclei) has enough proper time to reach the point under consideration. Since our approach is 
generally applicable for early (proper) times we can not unfortunately test for thermalization. In 
the best case one may follow the steps outlined in subsection [7.4| , conclusion 2., which provide a 
necessary but not sufficient condition whether a thermal medium is produced. The method, may 
also give an estimation of the times where thermalization occurs (if it does). We also illustrate 
how the momentum anisotropy and the spatial eccentricity may be obtained from the expression 
of T^ u (equation Q3.18D ) we have derived applying AdS/CFT and we compare our result with the 
results of the literature obtained by using pQCD or by hydrodynamical methods. In addition we 
argue about the necessity of an impact parameter and also about the fact that the bulk sources do 
contribute to the metric even at space-time points away from them. Furthermore, motivated by our 
result for the SE tensor, we propose how the energy density should evolve in more realistic cases. 
We construct a simple phenomenological model taking into account the Woods-Saxon profile of a 
nucleus at rest. We conjecture that the energy density at early times and close to the center of the 
collision should behave as in ( |7.29| ). Finally, we discuss the problem of negative energy density and 
how it is associated with the eikonal approximation in collisions. 



2. Setting up the problem 

We suppose that we have a single nucleus moving with the speed of light and use the AdS/CFT 
correspondence in order to model this in terms of a dual metric. We assume that the nucleus moves 
along decreasing x 3 and has a transverse (SE tensor) profile possessing azimuthal symmetry. The 
metric that realizes this situation corresponds to a metric 6 and a simple choice 7 may be taken to 

4 Certainly causality is already built to the dynamics of the problem. 

5 Totally there are six different areas but only three of them are not trivial. The other three may be obtained 
from the non-trivial regions by using mirror symmetry. 

6 As we see the nucleus is dual to a in gravity; hence we will use both notions interchangeably. 
7 In section || we will show how to construct more complicated metrics. These correspond to different nuclei 
profiles. 
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be 8 

L 2 

ds 2 = gMNdx M dx N = — \—2dx + dx~+ti(x + ,x 1 ,x 2 )z A dx +2 + dx 2 L + dz 2 },x M = x' 1 ,z. (2.1) 

z 

Here dx\ = (dx 1 ) 2 + (dx 2 ) 2 is the transverse metric and x ± = (t±x 3 )/\/2 where x 3 is the collision 
axis and t the time axis. With x M we include the four Minkowski components x^ = x 1 ,^ 2 ,^ 3 and 
x° = t and the component z which is the fifth dimension of AdSs with L being its radius. Now, t\ 
of fl2.1p cannot be an arbitrary function but should be constrained by Einstein's equations in AdS 
space. In the presence of matter these equations may be cast in the convenient form 

Rmn + J2 9mn = K 5 ( Jmn — - 9mn J J (2.2) 

with 

J = J™ = Jmn 9 MN 4 = 8nG 5 . (2.3) 
In equation (|2.2p, Rmn is the Ricci tensor, cjmn is the metric tensor in five dimensions and 



Jmn the bulk SE tensor 9 . G§ is the Newton's constant in five dimensions with [Gs]=-3. One can 
check that (|2.1| ) solves fl2.2|) exactly with a vanishing bulk tensor provided that 

' r3 d z -d 2 z -Vl)z%(x + ,x\x 2 ) = 0, (2.4) 



2 \^ 

where is the 2 dimensional Laplacian in the transverse plane. This means that the following 
condition should be satisfied 

-\z 4 Vlt l (x + ,x 1 ,x 2 )=0. (2.5) 

We are interested in rotationally invariant eigenfunctions for this differential operator. However, 
the only choice 10 is the \og(x 2 + x%) function which induces a delta function on the right-hand side. 
We interpret the appearance of such a delta function as a (nonzero transversely) localized Jmn- In 
particular, we choose the profile to be 

t 1 (x + ,x 1 ,x 2 ) = -fi\og(kr)5(x_), r = a/ (a^ 1 ) 2 + (x 2 ) 2 . (2.6) 



8 We choose the signature of the metric to be (— , +, +, +, +). 

9 The bulk SE tensor should not to be confused with the T M „ tensor of the produced matter in the four dimensional 
gauge theory. 

10 Apart from a trivial constant (see [92-94]). 
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In this equation k serves as a cutoff and its precise meaning will become clear in the next section 
while \i has mass dimensions [/i]=3. As we will see shortly, fi is associated with the (expectation 
value of the) SE tensor of the corresponding nucleus in the gauge theory and will serve as an 
expansion parameter of the metric. Taking into account that 

VllogO) = 2vr<5(f) r = (x\x 2 ), (2.7) 

we realize that there is a nonzero right-hand side for equations (|2.2|) associated with a five dimen- 
sional bulk SE tensor. It turns out that this bulk tensor is given by 11 

jW = ^ 4 dix 1 - b)6(x 2 )5(x + ). (2.8) 

All other components are zero. The presence of the superscript (W) on J ++ is to highlight that 
it is of first order in the parameter /i 12 . In this special case of a single , the first order solution 
happens to be the exact solution to all orders. In a diagrammatical approach, this solution would 
correspond to the diagram of figure [I]. It represents the measurement of the gravitational field at 
point x M , which, loosely speaking, is created by a single graviton emission from the bulk source of 
equation (|2.8j ) with effective coupling proportional to filog(kr). 

This bulk source is point-like along x 1 , x 2 and x 3 , has infinite extent along the z direction and 
moves towards decreasing x 3 . It represents a one-dimensional "string" of matter moving along the 
It is important to note that this bulk tensor is (covariantly) conserved as it should be. 
On the other hand, the boundary physics is obtained using holographic renormalization [136]. We 
find that the SE tensor corresponding to the metric (|2.1|) has only one non-zero component. To see 
this, we use 

where we assume that the metric is written in the Fefferman-Graham coordinate system 

ds 2 = ^ j^x", z)dx^dx v + dz 2 j. (2.10) 
Combining previous two equations we deduce that 

T 3 n (r^ 7~\ N 2 N 2 

(T++M> = T-PT lim / ++ V* J = -±t x = -^\og{kr)5{x + ). (2.11) 

n In equation ( |2.8| ) we assume that the center of the is shifted away from the origin by distance b along the positive 
x 1 axis. 

12 The superscripts emphasize the number of times the source ti (i=l,2 assuming we have two sources) appears. 
Generally in the object there exists the product t\t^~ k or any linear combination of differentiation/integration 
of t\ and t2 with respect to their arguments x^. 
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In both equations ( |2.9|) and Q2.ll ), N c is the number of colors. We note that T ++ (x >1 ) is positive 
for kr < 1 and therefore at first sight 1/k gives the maximum distance where the SE tensor of the 
nuclear matter is physical. In fact, this logarithmic (transverse) profile of the nucleus is a part of a 



more complicated profile which is both positive definite for all r and reduces to the one of ( 2.11 ) for 
small distances. Moreover, the corresponding metric is minimal in the sense of being the simplest 
one that can simultaneously capture some of the transverse dynamics of heavy ion collisions and 
allow for an analytical approach to the problem. We will give a detailed justification of these claims 
in the next section. 

Having defined all the necessary ingredients we now proceed to 
the main part of the setup. We want to superimpose two such shock 
waves whose sources are two one dimensional distributions of matter 
in the bulk. We want to collide them at a non-zero impact parameter 
and hence study the problem within the classical theory of gravity. 
Therefore, Jmn has in addition to (ET 



- /ilog(/.~r,;). i = 1,2 



the symmetric part 



J 



(i) 



7T/i 4 _ 
— jZ 0{X 



/->•- 



Six 1 + b)5{x 2 



(2.12) 



which creates a second . In terms of the gauge theory, this would 
correspond to colliding two nuclei in an off center process. Figure ^| 
represents the four dimensional picture, right before the collision of the 
two nuclei. Following [92], [93] and working with the Feffer man- Graham 
coordinate system, the metric that describes the process should look like 



ds< 



z 2 



2 dx + dx + dx\ + dz* 



(x H 



b, x 



Figure 1: The solution: a 
graviton is emitted from the 
bulk source with coupling 
/ilog(/cr) which is measured 
at the point x . This is 
a very special case where 
a single graviton exchange 
between the source and the 
bulk happens to be an ex- 
act solution to the nonlinear 
Einstein's equations. 
z A dx +2 



+ t ( i\x~,x 1 + b,x 2 )z 4 dx~ 2 + e{x + )6{x-)g {2 }(x K ,z)dx' i dx u + ... 
t^l(x l Tb,x 2 ) = -/ilog (k ^(x 1 T b) 2 + (x 2 ) 2 ) 6(0*). 



(2.13) 



The first three terms correspond to the empty AdSs. The next two are of first order in /i and are 
created by the two lines of distribution of bulk matter that stretch along the z direction. These 
objects move (initially) towards each other along x 3 and they have an impact parameter 2b along 
the X cLXIS clS figure |] depicts. As they are the sources of the two shock waves, they correspond 
to two vertex diagrams that look like the one in figure [l]. This is a superposition of two metrics 
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~ /ilog(fcri) 



Figure 3: The Feynman diagram that represents the fi 2 correction of the metric: It represents the 
first nontrivial correction to (2.13). It shows how the two metrics each of which like (2.1) merge. The 
gravitational field is measured at the point x M . 



with each one looking like ( [2.1| ) while the gauge theory picture is shown in figure |2|. However, the 
non-linearities of the gravitational field require higher order terms. The second order corrections 
are explicitly displayed in fl2.13| ) and they appear once the two shock waves cross each other in 
the forward light cone. This is precisely the meaning of the theta functions; they emphasize that 
the metric Q2.13 ) solves Einstein's equations exactly in the presence of both shock waves only for 



negative x . According to holographic renormalization [136] and equation ( 2.111) , the additional 



terms of the metric that appear in the forward light cone correspond to the matter of the produced 
medium in gauge theory after the nuclear collision. The main work of this paper is to calculate them 
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to order /z 2 , that is find gffi and determine of the gauge theory. The second order correction in 
\i of g^ u corresponds to the diagram of figure ||. 

3. Choosing the transverse profiles of the initial nuclear matter 

In this section we justify the choice of the profiles of the shock waves which look like t\ of equation 
(f2.6p. According to equation and ( [2.7| ) we have that 

- z & z - d 2 - ViJ {-z 4 log(Ax)) = 2irz 4 5{r) (3.1) 

We will see this differential equation from a different perspective. We will search for the Green's 
function of the operator ( |2.4|) and then use it in order to determine t± by assuming that the source 
of this operator behaves as ~ z 4 5(f) . We need the function G which satisfies 

fh z - d 2 z -V 2 ^jG = 5{z - z>)5{r- f ) (3.2) 

and which is finite as z — > 0, oo. We will see that the solution ti ~ z 4 \og(kr) is only a part of a 
more complicated metric and we will also associate k with an IR cutoff of the fifth AdS coordinate 
z. Searching for solutions of the form 

G = J ^ 2 e^') A (q,z } z') (3.3) 



we find that A should satisfy 



(l dM - Ql + q ^ A = 5(z - z') 



(3.4) 



It turns out that the two independent solutions for A are z 2 K2(qz) and z 2 l2(qz) where K2 and I2 
are Bessel functions. Now, the boundary conditions together with continuity of G require that A 
should behave like 



A ~ 6{z' - z)z 2 K 2 {qz')l2{qz) + 6{z - z')z 2 K 2 (qz)I 2 (qz'). (3.5) 
In order to determine the overall coefficient, we substitute (|3.5|) in (|3.4|). The result is 
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z 2 



A= — [6{z' - z)K 2 (qz')I 2 (qz) + 9{z - z')K 2 (qz)I 2 (qz')} . (3.6) 
The next step in to plug ( |3.6|) in (|3.3| ) and perform the angular integration. This yields to 



2 i-oo 

G = - — - / dqqj (q\f- ?\) (6(z' - z)K 2 {qz')I 2 {qz) + 9(z - z')K 2 (qz)I 2 (qz')) . (3.7) 
Anz Jq 

These integrals 13 are tabulated in [171] and the answer involves the associated Legendre functions 
of the second kind (Q„). We find 



* e^ 2 -l)-iQ^ H; u= \r-rf + (z-zr + ^ (g 



G 



The function ma y be simplified even further [171] in terms of a hypergeometric function. The 
final result is 14 

G(f, z- P, z>) = y^^F(2, 3/2, 3, l/u\ (3.9) 
We may also cast last equation in the rather simple form 



if -> / 2z 

G(r,z;P,z') = — — - 4|f - P\ 2 + 4z 2 1 



/2 



f - r'| 4 + (^ 2 - ^ /2 ) 2 + 2r 2 (z 2 + z' 2 ) 



+ 4 \^-z' 2 + \J \r - r'\ 4 + (z 2 - z' 2 ) 2 + 2r 2 (z 2 + z> 2 ) j j. (3.10) 

This result agrees with [114], [115] and [142] and it was obtained more elegantly by using the 
symmetries of the AdS^ space and reducing the problem to solving a differential equation of a single 
variable u. 

Now we will make use of the Green's function we derived in order to answer the question we 
posed at the beginning of this section. We wanted to find the solution of the following differential 
equation 



13 Both integrals converge as q — > oo because of the presence of the 8 functions. The values of the integrals turn 
out to be the same under z O z' and hence independent of the ordering of z and z' . 



The fact that G of (3/7) is not symmetric under z <-> z' is because the differential operator ( 3^2 ) is not Hermitian 
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d z - d 2 - Vi )h = 2nz 4 5(r- P)5{x+). 



(3.11) 



We have to take the product of the right-hand side of (|3.11| ) with ( p.lO|) and integrate. The 
transverse integration is trivial. But then, one may observe that the integrand of the product of 
( P-llj ) with (|3.10|) at large z' behaves as Zj. The only available scale that fixes the dimensions is r 
and hence the result should receive a z 4 log(— ) contribution from large z' . This suggests we should 
place an infrared cutoff z c on z' . 

Equivalently, we may argue in a more physical way. We may alternatively ask for the solution 

of 



z 



d 2 



h = 2ttz 4 8{z c -z)8(r- r')5{x + ) 



(3.12) 



assuming that the bulk SE tensor does not extend to infinity in the z direction but the source in 
the bulk looks like a rod of finite length z c . The exact result for h for this bulk tensor is then 



h = G ® 2nz' 4 9(z c - z')S(f- r)5(x + ) 




2(r 2 + z 2 )z' 2 - z' 4 + (r 2 + 3z 2 + z' 2 )^^ + (z 2 - z' 2 ) 2 + 2r 2 (z 2 + z 



/2^ 



+ 4z 4 log (2 (r 2 -z 2 + z' 2 + Vr 4 + (z 2 



z' 2 ) 2 + 2r 2 (z 2 + z' 2 




-^)(log(^) + ^-^) + 0(l), 



(3.13) 



up to the term — |z 4 5(x + ) which we have absorbed in the cutoff z c . Our approach to the problem is 
to work only with the logarithmic piece of ( |3.13| ) assuming that z c is large enough. Thus, whatever 
results we will obtain using the shock waves of ( [2. 13] ), should be thought of as arising from the 
(exact expression of the) shock waves of ( |3.13| ), up to corrections of 0{\). We stress that the 
logarithmic piece is also an exact solution of Einstein's equations with a bulk stress energy tensor 
that has the form of ( |2.8| ). This bulk tensor is different than the right-hand side of equation ( 3TT2" ) 
but the two agree in the limit where z c — > 00. 

We thus observe in a direct way how the z direction of AdSs plays a crucial role in the transverse 
behavior of the gauge theory SE tensor [154] - [170]. In addition, as r — > z+ the gauge theory SE 
tensor of the that corresponds to the exact expression for h, equation ( |3.13| ), remains positive. In 
order to see this, we expand the exact solution of (|3.13|) to order z 4 obtaining 
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I-, ^ / Is _ 2x 2 + 3 \ 4 r 



6(x + ) 41og(l + -)-2-— — - \z\ x = -. (3.14) 
\ x z (1 + x z ) z / z c 

Equation (|3.14j) is (proportional to) the exact T^ u of the gauge theory which corresponds to the 
exact produced by a rod of length z c according to the right-hand side of ( 3.12|) . It is monotonically 



decreasing to zero and is strictly positive for all r. It reduces to the approximate expression of 
( |3.13| ) when z c — > oo. It is now evident that we may identify l/z c with ~ k of equations fl2.6|) and 

(El)- 

In addition to the above analysis, we may derive a similar conclusion by working from a different 



perspective. We may associate the (IR) cutoff of the parameter q (see below) of ( |3.7| ) with some 
infrared scale, probably comparable with Aqcd- We therefore convolute G of ( |3.7| ) with the right- 
hand side of ( [3. lip and perform the trivial d 2 r' integration and also the z' integration. We obtain 

h = G® 2nz' 4 6(r- 7)5(x + ) = z 4 [°° dq J °^ T - ^ 6(x + ). (3.15) 

Jo Q 

We immediately encounter a singularity at q — 0. In order to regulate the singularity, we place an 
IR cutoff on q. Now, by performing the q integration we obtain the final result 

ti = = / d <i s ( x ) 

z Jk H 

( 1 2 2 / /22\ ^ 

= | - log(A:r) + —^—FppQ (l, 1; 2, 2, 2; + const. | 6(x+) (3.16) 

where Fppq is a generalized hypergeometric function. Therefore, as long as we work in the kine- 
matical region r, r < 1/fc we may ignore the hypergeometric piece of equation ( |3.16| ) and safely 
work with the logarithmic piece 15 . 

Working with the logarithmic piece of the transverse profile means that the energy density of 
the single beam moving along the light cone is infinite but integrable at the center of the beam and 
tends to zero logarithmically and radially outwards; it vanishes at r = 1/k. This gives an additional 
indication that 1/k should be of the order of the nucleus size [154] - [170], that is 

k ~ ^ (3.17) 

As 

where A is the atomic number. Under this assumption, the total energy of the beam can be 
estimated by 



15 See section ^ for details. 



13 



V* . _ A 



E ~ — /x / dx + / dr5(x + )r log(fcr) ~ /x— ^ . (3.18) 

</-oo JO ^-QCD 

Taking into account that the energy of the matter moving along x~ should be proportional to the 
p" momentum of each nucleon of the nucleus and hence E ~ p~A, we deduce that 



V~p-A 2 Q CD^ (3.19) 

which is in agreement with [92,93]. What is more is that this choice of a nuclear profile is very 
simple in the sense that the z and the transverse coordinate dependence r of the factor. This will 
simplify the calculations significantly although they are still not trivial. We note that we were able 
to construct nuclear profiles that behave like ~ — with n being an integer. These correspond to dual 
metrics that are much more involved and the calculations would have been harder. Finally we find 
it attractive to work with a metric that replaces the ~ — z A \og(kr) profile with ( |3.9D . This metric 
gives a nuclear profile 

^ ( 2 t" /2\3 [H^l] with being an arbitrary constant. The advantage of 

such a metric is that it gives a point-like bulk (five dimensional) SE tensor and the nuclear profile is 
everywhere positive. Unfortunately this metric is relatively involved and makes calculations hard. 
We leave the investigation of such a case for a future project. 



4. Back-Reactions 

4.1 Corrections to Jmn and Geodesies 

As has been already mentioned below (|2.8| ), is conserved in the gravitational field of Q2.1|) . In 
fact, conservation for this case happens to be valid to all orders in /x 16 

V M jfiJv = 0. (4.1) 

Conservation to first order is still valid when we consider simultaneously the bulk distributions 
of matter J+ 1 } and J__ in the presence of the gravitational field fl2.13|) . However, this is no longer 

true at the second order in /x. The reason is because the J ( J++) source moves in the gravitational 

field of the t\ (i 2 ) , altering its initial trajectory. Figure (| outlines what happens while figure |5] 
offers a diagrammatical intuition regarding the self-corrections to Jmn- This implies that we should 
correct J MN in order to preserve conservation (of the bulk Stress- Energy tensor). However, since 
we do not know the nature (equation of state) of Jmn we make the assumption that these objects 
interact only via gravitational forces 17 . 



16 In practice, only the first order in /i appears in the resulting equations. This is in accordance with our intuitive 
picture of figure El Gravity behaves linearly with respect to the metric (2.1). 



17 



More precisely, we assume that any other interactions are small compared to the gravitational forces. 
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Figure 4: The sources (one dimensional distributions of matter) represented as vertical lines extending 
along z and their trajectories (suppressing the x 2 -axis): For negative times they move along straight 
(dotted) lines. At t=0 each distribution intersects the shock due to the other distribution and its trajectory 
suffers a kick (see curved path); hence Jmn changes with time. 



J++ J++ 

„M 



' a ' l M/w\J» 



x M 



Figure 5: Graviton emission diagrams resulting from self corrections to Jmn- The first source interacts 
via the gravitational field created by the other and vice versa. The point x M is the space-time point where 
gMN is measured. 



4.2 Calculating the Corrections for Jmn 

If the bulk matter we have to deal with consisted of just two point-like particles, they should have 
traveled along geodesies; as it is rigorously shown in [172] conservation of (the total) Jmn is then 
guaranteed (for an application see [175]). In this case we would need 
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J MN = mJ2 *(/)^)7^ (4) (%)-%)(«(/))) * M = {x\x\x*,x\x\z) (4.2) 

(7=1) V 

where m has dimensions of mass while x includes the four spatial components. Equation ( |4.2j) gives 
the total bulk SE tensor of both point-like particles each one moving along the trajectory xm(sm) 
parameterized by sm. The quantity g is the determinant of the (total) metric tensor while the 
dots denote differentiation with respect to the parameter sj. However, in our case we deal with one 
dimensional distributions of matter extending along the z direction. Hence, we naturally generalize 

(op to 



jMN = ^y" xf/ } x^-=5 (3) (x (I) - %)(s ( j))) p MN (z) (no summing over M,N) (4.3) 
«5 (/=1) V 9 

where the factor ^ has dimensions of mass, is first order in fi and its exact form is a convenient 



5 (/=i) 

convention (see below) while x^ does not include the z coordinate 18 . In this equation we assume 
that the distribution of matter is made from point-like particles distributed in a continuous way in 
space-time along the z-axis. The total bulk SE tensor is the summation (integral) of all of these 
particles which are assumed to be glued tight together so that there are no mutual interactions. 

Under these assumptions, one may compute the 0(fi 2 ) corrections to J MN from both distri- 
butions. Before doing so we find it instructive to check whether this formula reproduces (|2.8|) in 
the case of one distribution (/ = 1) which moves ultra-relativistically along negative x 3 and is 
parameterized by 



x 



(\\{x ,x ,x ,x ,^) — (t, — /■, b. (). z) . 



(4.4) 



In light cone coordinates and choosing to parameterize the trajectory by x , equation ( |4.4| ) then 
implies 



x (i) = (x + (x ) = 0,x (x ),x 1 (x ) = b,x 2 (x ) = 0, z(x )) = (0,x , b, 0, z). 



Direct substitution of ( fig ) to ( fOQ yields 

J (i) = J (i)++ 



^^ 4 5(x + )5(x 1 - b)5{x 2 ) 



(4.5) 



(4.6) 



for p (z) = J dz' ' yj —g(z')z' i 5(z' — z) = ^—gz A reproducing equation ( |2.8| ). This computation also 
clarifies the convention of ^ in (|4.3|) . The second piece of J^, equation ( |2.12|) , due to the second 
particle is reproduced similarly. 



18 In contrast with (E 
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The next step is to calculate the next (second) order corrections (in fx) of J^ u which translates 
into finding the corrections to the trajectories xm. These may be obtained from the geodesic 
equations. In particular, we only need the first order corrections to xm as we already have a power 
of \i in front of the summation operator (see ( |4.3| )). The geodesic equations we need read 



x M i + ry. NP xfxf = I, J =1,2 I^J (4.7) 

and are interpreted as the motion of distribution I in the gravitational field of the distribution J 
(due to Tj. NP where V are the Christoffel symbols) and vice versa; this is precisely the meaning of 
the subscripts I and J. We begin with computing the corrections to the distribution 1 = 1 whose 
(first order) perturbed trajectory looks like 19 

(x M (x-))(°) + (x A/ (x"))« = ({x+Y 1 \x- + (ar)« b + (x 1 )^ (x 2 ) (1 \ (4-8) 

where we have chosen to parameterize the trajectory with x~ (i.e. S(/=i) = x~). The superscript W 
denotes the order (in fj,) in the expansion. Taking into account (|4.3|) and the fact that T^ NP ~ A* 
otherwise is zero, we deduce that both of the terms x N or x p have to be of zeroth order; i.e. the 
only choice is iV = P = —. This implies that we need to determine T__ to first order in /x that 
arise from the second particle 20 and which read 

4 4 4 

rl_ = -yV- r:_ = o v l __ = --t 2 ^ r 2 _ = --t 2>x2 r*__ = -z%. (4.9) 

A few explanations about our notation are in order: the term t 2 is due to the second particle and is 
given in equation ( |2.13| ). The subscript 2 ^ on t 2iXl ± denotes ordinary differentiation of the source 21 



t 2 with respect to the coordinate x M . According to (|2.13 ), the source t 2 is of first order in /j 22 and 



as a result, the same applies for the T's of (|4.9| ); it should be by now obvious that the Christoffel 
symbols are due to the second distribution (i 2 ). 

The final step is to integrate ( |4.7|) using (|4.9|) and using causal boundary conditions and substi- 
tute the result in formula ( |4.3|) . As we are interested in second order corrections in /i, we immediately 
conclude that at least one of x M or x N should be of order zero, i.e. M = — , N ^ — or M ^ - , 
N = — . We also note that \J—g ~ 0(1) + 0(/i 2 ) and hence according to ( |4.3| ) corrections from 
\J—g do not contribute to Jmn at 0(/i 2 ). We consider two cases 



Case I: J=g~ = 0(1), M = N = - 

19 We drop the subscript m which labels the particle 1=1 for simplicity. 
20 Where we (also) dropped the subscript / 2 ) which labels the second particle. 

21 As matter implies curvature and curvature implies matter, we will often use these notions interchangeably while 
their distinction should be evident from the context. 

22 From now on we drop the superscript ^ which denotes the order in ^ from t^\. 
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In this case the modification of Jmn of equation (|4.3|) to the order we are working is in the 
arguments of the delta's. We have 23 



J =—^z A 5 ( x + / dx~ 



z 
~2 



dx t 



2,x~ 



5 { x 1 — b / dx 

2 



dx t 



2.x 1 



x 5 \^x 2 — — J dx J dx t 2)X 2 J + the integrated 8{z) term (4-10) 

where by "the integrated 5(z) term" we mean the contribution due to the deflection along z in- 
tegrated twice: the first integration takes into account that the point particle of the distribu- 
tion 1 = 1 located at z, interacts with the matter distribution J = 2 and then the second 
integration sums over all particles of "type" 1 = 1. It turns out that this term is equal to 
^z 6 J dx~ J dx~t 2 8(x + )5(x 1 — b)5(x 2 ). Plugging this term in ( (4.10|) and expanding the delta's 
to first order in the sources we obtain 



J- 



7T/X 4 



- b)5(x 2 ) dx \ t 2 5'(x+)5(x 1 - b)5(x 2 ) 



+ I dx- (t 2 ^5(x + )8\x l - b)5(x 2 ) + t 2 ^5(x + )5(a; 1 - b)5' (a; 2 )) 



- 8z 2 I dx~ I dx~t 2 5{x + )5{x l - b)8(x 2 ) 



(4.11) 



We may cast (the first order correction terms of the) last equation in a compact form by expressing 
it in terms of t\ and t 2 . Using the identity 



/x5(x ± )5 (2) (f-6i j 



2tt 



Vi*i, 2 6 li2 = (±6,0) 



(4.12) 



(see (O) and (513)) the 0(/i 2 ) terms of (KW) take the form 



(^(D ++ ) (2) = (W (2) = ^72 I dx~ 



x (^z 8 (t 2 V 2 ± t 1;X + + V 2 ± t 1>x i J dx~t 2<x i + V 2 ± t hx 2 J dx~t 2jX ^j + 16^ 6 V 2 L ti J dx~t 2 j (4.13) 

where we restored the subscript (i) and the superscript ^ in order to highlight that this is the 
second order correction to Jmn of the first distribution. 



23 All the integrations with respect to da;* that follow from now on (see also Appendix |a|) will imply the obvious: 
f dx~ stands for f x dx'~ , f dx~ f dx~ stands for f x dx'~ f dx"~ etc. 

J J—oc ' J J J—oo J— oo 
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Case II: y/=g = 0(1), M = N ^ - 



In this case the modification of Jmn of equation ( [4.3[ ) to this order we are working is in the 
factor x M x~ = x . Combining Q4.7|) and ( |4-9|) one may compute x . Plugging this result to 
and employing the identity ( [4. 12| ) in order to write the transverse delta's in terms of t\ yields 



(4.14a) 
(4.14b) 
(4.14c) 
(4.14d) 



( J JD (2) = - 


V 2 ^ 1 


(</(i )+ i) (2) = 


-(j- 1 )^) = 


(^(D +2 ) (2) = 


-(rf = 


(</(d+,) (2) = 


_(J-*)(2) = 



4k§ 



rVlti / dx~t 



Second order corrections to the total J, 



MN 



The second order corrections to the stress energy tensor ( J(i)mn)^ of the first distribution is 
given by equations (|4.13|) and ( [4.141) . The corrections (J( 2 )mn) of the second distribution may be 
found analogously and therefore the second order corrections to the total Jmn read 



(J + -) {2) = -(J + -) {2) = (hV 2 ±tl + tl vlt 2 ) 

4 K 5 



(Jn) {2) = (J 22 ) {2) = (Ji 2 ) {2) = (J zz ) i2) = o 



(4.15a) 



x \^z 8 \t 2 V 2 ± t 1>x + + V 2 J ltX i J dx~t 2>x i + V 2 Ji )X 2 J dx~t 2tX 2j +16/Viti J dx~t 2 J (4.15b) 

x (* (t^- + Vi^ + V|* 2 , J&tJ) + 1^/^) (4,5c) 

( J +i) (2) = 7^2 V i^ / ^a,* 1 (J-i) (2) = ]\^lh [ dx + t 1>x i (4.15d) 

(J +2 ) (2) = j^vit! / dx~t 2 , X 2 (j_ 2 ) (2) = -^v 2 ^ / dx + t liX 2 

{J +Z )V = ^Vit! / dx-t 2 {J^ = ^Vit, / dx + h 



(4.15e) 

(4.15f) 
(4.15g) 
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The first equality of (|P5a|) is not completely obvious and so we prove it below by considering 



(J + -) {2) = (g+ M g- N J MN ) m 



2k\ 



viti ) + {z%){-\) 



24 



vl*i 



(-i)(-i)^ 8 (-^(^vi^ + ^vito) 



1 1 



-— z 8 (hVlh + hVlh) 



(4.16) 



where in the fourth equality we used the fact that the (bulk) SE tensor of the point particles (see 
( p.8| ), ( |2.12|) and ( |4.12| )) to first order in \x takes the form 



J 



(i) 



24 



24 



(4.17) 



Although we will mostly work with (|4.15 ) which is a compact expression, for concreteness, we write 
Jmn of ( |4.15| ) in terms of the coordinates in order to clarify its form and show explicitly that it is 
localized (in all but the z directions). Defining 



f[ = f — b\ f*2 = r — 6 2 - 
and employing ( |4.12| ) in ( [4. 15| ) we obtain 



(4.18) 



(j + _)( 2 ) = ^z s log(2kb)5(x + )5(x-) (5 (2) (n) + <5 (2) (r1))) 



24 



,(2) _ TH£ 



24 



z 8 log{2kb)5'{x + )5 i2 \r 1 ) + z 8 -^—^5{x + )5'{x 1 - b)5{x 2 ) 



_ 2 

h/ , TO X . % ^J (x + )S(x l -b)5'(x 2 ) + l6x-z 6 log(2fc|6|)5(x + )5 (2) (r"l) 



7T/i 

Anlb 



Ab 2 + (x 2 ) 2 
2 



z 8 6(x-)5(x + )5 {2 \r-i), 



(J + 2) {2) = 0, 



(J- 



»(2) 



2n/j, 2 7 



z 7 d{x-)\og{2k\b\)5{x + )5 {2 \r l ). 



(4.19a) 



(4.19b) 

(4.19c) 
(4.19d) 
(4.19e) 
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The asymmetry between J±l and J+\ is due to the fact that the impact parameter b has only x 1 



component (see (|2.13| )). The remaining nonzero components that complete (|4.19|) may be obtained 



using the symmetries of the problem: J^_ and may be obtained from J+\_ and respectively 

(2) 

by interchanging + -h- — and b -H- —b. We observe that J MN is localized in all but the z direction 
and exists for positive times or exactly at the collision point at the collision time (this is the j| 2 2 
component-see ( |4.19a| )). 

4.3 Conservation, Tracelessness and Field Equations 

The second order corrections to (the total) Jmn have already been calculated in the previous section. 
One may check that the action of a covariant derivative (contracted) on Jmn yields 

( (V M)(o) + (V M)(D) ((J MJV )« + (J MN )W) ~ 5 N± VlhVlt 2 + O(^) (4.20) 

where V denotes a covariant derivative, 5^± is the Kronecker delta while the superscripts denote 
the order in \i. According to (|4.12|) the expression V^tiV 2 ^ gives 

VlhVlh ~ 5(f- Si)«J(r - b 2 ) ~ 5(b 2 - b x ). (4.21) 

Therefore, conservation of Jm± (to 0(fi 2 ) at least) is valid if and only if the impact parameter is 
not zero. We thus see that introducing a nonzero impact parameter is essential in order to preserve 
conservation. This is one of the important conclusions that we derive in this work 24 . 



It is also useful to compute the trace of Jmn as it enters the field equations (see (|2.2| )). A short 
computation yields 

J = 9 MN Jmn = ( 9 MN ) {1 \Jmn)W + (g MN r\J M N) (2) = + 0(/i 3 ) (4.22) 

which shows that the stress-energy tensor is traceless to order [i 2 . Tracelessness is very convenient 
as it simplifies Einstein's equations which become 

Rmn = 4Jmn + 0{^) 4 = 8ttG 4 . (4.23) 

5. Field Equations 

In this section we will write down the field equations ( |2.2| ) to order fi 2 . These determine the 
functions gffi 25 of ( |2.13| ) which is what we need to calculate. We will write down their linearized 



form expanding around the background that is determined by the simultaneous presence of both 
shock waves and show how they may be decoupled. We will find out that g^J obeys the equation 



24 Later we will see that a as b — > we have to deal with another problem: the metric tensor also diverges. 
2b g$ are not exactly a subset of the components of gMN since according to their defining equation (2.13), they 

2 

are off by a factor of 4^ . 
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= (5.1) 

where □ = (-d z — d z — + 2d x +d x -) is the D'Alambertian in AdSs and s^ v some tensor related 
to J^y of equation fl4.15| ). Although decoupling all of the metric components is possible and has 
been done, we only display the components of ( f4.23p that we will use. We need 26 



v++) 7r[3g++,z - zg++, zz - zg ++x 2 x 2 - zg ++jX i x i 
2z 

+ zg+2,x+x 2 + zg + i tX + x i — zgn tX + x + — zg 2 2, x +x+] = k\.J ++1 (5.2a) 
[-\ — ) — [ — lQz z tit 2 + z (t lx +t 2x - — 2U x it 2x i — 2t lx 2t 2x 2) 

- 2g n>z - 2g 22tZ + Wg+~ z - 2zg + _ :ZZ - 2zg + _ iX i x i - 2zg + _ :X 2 x 2 

,x~x z + zg+i :X - x i + zg_ 2 + zg 

- 2zg ++>x - x - - 2zg__ :X+x + + Azg + _ :X+X - - 2zg u>x + x - - 2zg 22 , x + x -) = k 2 5 J + -, (5.2b) 

(+1) -r[z 9 t 2 . x it l x + + 6g + i, z - 2zg +lzz - 2zg +lx 2 x 2 + 2zg +2 x i x 2 - 2zg ++ x - x i 
4z 

+ 2zg 12jX+x 2 + 2zg + _ )X+x i - 2zg 22}X + x i + 2zg +1}X + x - - 2zg_ 1}X + x +] = «f J+i, (5.2c) 

(11) 2~[ - 8z 7 t 1 t 2 + z 9 (t 1>x it 2)X i + t ltX i x it 2 + ht 2iX i x i) 

- '2g+-,z + ^gii,z + g22,z - zg u>zz - zg xltX 2 x 2 + 2zg 12 ^ x i x 2 

+ 2g + _ x i x i - zg 22>x i x i - 2zg +ltX - x i - 2zg^ x+x i + 2zg 11>x + x -) = k\j xl = 0, (5. 2d) 

(12) (t 2iX 2t 1>x i + t 1)X 2t 2)X i + 2t 2 t ljX i x 2 + 2t 1 t 2jX i x 2) + 6g 12jZ - 2zg 12tZZ + Azg + _ :X i x 2 

- 2zg +liX - x 2 - 2zg +2>x - x i - 2zg_ hx+x 2 - 2zg_ 2)X+x i + &zg 12iX+x -} = K 2 5 Ji2 = 0, (5.2e) 
(zz) ^-[32z 7 tit 2 - 2# + _, 2 + g n>z + g 22)Z + 2zg+- tZZ - zg lljZZ - zg 22)ZZ ] = k\.J zz = 0, (5.2f) 

(lz) ^[Qz 7 t hx it 2 + 6z 7 t 2>x iti + g 12>x 2 z 

+ 2g+-,x i z - gm^z - g+i,x- z - g-i,+ z ] = k\J Xz , = o (5.2g) 

(2z) \- [6z 7 t 1/x 2t 2 + §z 7 t 2iX 2t x + g 12tX i z 

%g+-,x 2 z - gn,x 2 z - g+2,x- z - g-2,x+ z ] = k^J 2z = o, (5.2h) 

+ ,x~z 

+ g+-,x+z — gu,x+z — g22,x+z] = k\J+ z (5-2i) 



[+ z ) ^z 7 t 2 t ltX + + g +2 ,x 2 z + g+i^z - g++, 



26 From now on, we will drop the superscript (2) (which denotes the order in the expansion in fi) on gffi and jffi 
for simplicity. 
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where the components of Jmn are given by ( |4.15|) . In order to decouple these equations we have 
to take appropriate linear combinations of derivatives or integrals (with respect to the coordinates 
) of several of the components of ( |5.2| ). As a first step we integrate ( [5.2f| ) twice with respect to 



x 



M 



z from zero to z with the boundary condition that vanishes as z — > 0. In particular, we apply 
the operator z 2 /* 1/z' 3 /*' 2z"dz"dz' in both sides of {z"z") of ( ggjj ). The result is 



-2g + _ + g u + g 22 = -z 8 ht 2 (5.3) 

This equation implies that the z 4 dependence of this equation is on the left hand side only; accord- 
ing to ( |2.9|) this dependence is proportional to Tff. Therefore, equation ( |5\3| ) is the tracelessness 
condition of the T^ u gauge theory tensor. 

We now verify that conservation of is compatible with ( |5.2| ). Let us consider for instance 
( |5.2fJ ) and find the z A dependence of this equation which is given by 

(gi2,x 2 + 2flf + _^.i - g22, x i ~ 9+i,x- - 9-i,x+) \ z 4 = (5.4) 
and which together with the tracelessness condition (— 2g + _ + gu + #22 1 4 = 0) implies 

(0i2,*a + 9n,x^ - 9+i,x- - 9-i,x+) \ z a = 0. (5.5) 

But equation ( |5.5|) is exactly the conservation equation for T^ u at order /i 2 , that is 

(V%,) (2) =0 v=\. (5.6) 

In the same way one may see that (V^T^,,)^ = for each v = +, — , 1, 2 and conservation is encoded 
in the [yz) components of the field equations and they hold order by order in //; consequently we 
have 



ij,i> 



(5.7) 



to all orders. 

Next, we show how equations Q5.2 ) may be decoupled. For this, we consider the linear combina- 
tion (11)— 29^1 (Jj (lz')dz') 27 which involves the components (5.2d) and (|5.2g ) and employ equation 
Q5.3D in order to get rid of the linear combination 2g X2 + g 22 that appears on the way. We apply a 
similar method for the g 22 component. For the g\ 2 we take (12) — d x 2 (J* (lz')dz') —d x \ [J^(2z')dz ,s j 
where we make use of equations (|5.2ej) , ( 5.2g|) and (|5.2h|) . The is decoupled by considering the 



combination (H — ) — d x + (f*(—z')dz') — d x - (f Q (+z')dz') which makes use of (|5.2b|) , ( |5.2i|) and the 



27 With this notation we mean that we integrate and differentiate ( |5.2g ) with respect to z and x 1 respectively and 



subtract the result from equation (5. 2d). A similar meaning applies to all combinations that follow. 
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corresponding equation for (—z) 28 . For the g ++ component, we take (++) — 2d x + (f Q (+z')dz') 
where we use equations ( |5.2a|) and ( |5.2i|) . Finally for the component g + \ we make use of the equa- 
tions ( |5.2c| ), ( |5.2g| ) and ( |5.2i| ) by considering (+1) — d x + (f(lz)dz) — d x i ( K §(+z)dz) with a similar 
combination for the (+2) case. In all, cases, we make use of ( |5.3| ) in an appropriate way. After 
some algebra and using (|4.15|) we eventually obtain 29 



\-d z - d z - V ± + 2d x +d x -\g u = -z *i*2 + z [^t hx it 2jX i - -t hx i x it 2 - -tit 2<x i x i J (5.8a) 

0d 2 - d 2 z - Vi + 2d x +d x -^jg 12 = ^z 8 (2t 2>x 2t 1>x i + 2t 1>x 2t 2>x i - t 2 t 1>x i x 2 - ht^i^) (5.8b) 
0<9 2 - d 2 z - Vi + 2d x+ d x -^jg ++ = -^z%t ltX+x+ + 8z 6 Vlh J dx~ J dx~t 2 

+ ^z 8 J dx~ J dx~ (t 2>x iV 2 ± t 1>x i + t 2:X 2V 2 ± t 1>x 2) (5.8c) 



(5.8d) 



Q<9 2 - d\ - Vi + 2d x+ d x -Jg +1 = \*?(- h,x^ x , x+ + 2t 2 t liX+x > 

+ ^z 8 J dx~ (t^iVlh - t 2 V 2 ± t ltX i) 
(jz^ z ~ ^ ~ ~^ ^x+d x -jg+- = — —tit 2 z 6 + z 8 (^-ti ;X +t 2jX - — t ljX it 2)X i — t ltX 2t 2>x 2 

- \hVlt 2 - ^ 2 Vit x ) (5.8e) 

The remaining differential equations (and hence solutions ) of g^ v can be determined as follows: 
g^ 2 are obtained from by interchanging the coordinates 1 <H- 2 while the components g_^ are 
obtained from components g +At after interchanging simultaneously x + <H- x~ and bi <H- b 2 where b\ t2 
are given by ( |4.12| ). These equations will be solved in Appendices [A] and |B] 30 . 



Remark: 



The (z 4 coefficient of the) component g + _ may and will be obtained from Q5.3|) and so we 
will not solve (|5.8ej ). However, as a cross check we verify that ( |5.3|) which has been obtained 
from (|5.2fl) , provides the same information as (|5.8e|) . In order to see this, one may simply add 



equations ( |5.8a| ) and the corresponding differential equation for g 22 and subtract twice ( |5.8e[ ) yielding 
-d z — d 2 — Vi + 2d x +d x - I (gn + g 22 — 2g + _ — ^zHiti) = 0; this justifies the claim. In fact, this is a 



28 It may be obtained from (+z) using the obvious discrete symmetries of the problem. 

29 The corresponding solutions for the components g_ M , /i = —,1,2 are obtained from g +ll by interchanging the 
coordinates + f-> — and setting the impact parameter b — > —6. 

30 More precisely, we need just the z A coefficient of and we will derive only this part of the metric. 
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general feature of counting the degrees of freedom of both sides of the AdS / CFT duality correctly: 
there are 10 degrees of freedom for in the gauge theory side while in principle in the gravity 
side there 15 degrees of freedom for the metric that are specified by 15 differential equations that 
result from the field equations. Hence there should be a redundancy. Obviously the redundancy 
for the metric comes from the coordinate choice freedom while 5 out of the 15 components of the 
field equations should be satisfied trivially. Indeed, in the Fefferman- Graham system we are using 
here, the trivial components of the field equations (|4.23 ) are the (zM) components of ( |5.8|) and 



they encode, as already been said, the conservation and tracelessness of T, 



6. Deriving T^ v and results 

6.1 Deriving the Retarded Green's function 

We wish to calculate the components g^ u and of G5.8|) . In order to do this, we need to invert 
the differential operator (-d z — d 2 — + 2d x +d x -). In particular, we are looking for the Green's 
function of this operator with the boundary condition that it vanishes at z = and also at negative 
times t; that is we are looking for the retarded Green's function. This has already been done in [173] 
31 and independently but later in [92]. The integral representation of the retarded Green's function 
then is 



oo 

z 2 



G(t,x, z;t' ,x , z) =6{t — t')— I dm 



z' J J (2tt 



d 3 k 







sinUt - i!)\/m 2 + k 2 ] T , s T/ , N ,„ 

x m — ' -e lk[x x > J 2 {m z) J 2 {m z ) (6.1) 

ym 2 + k 2 

where 3 ). Performing the angular integration of the k 32 variable yields 



oo 

r ^ z 2 1 ^ I A _ f°° dk m Sin[(t-f)Vm 2 + k 2 } 



z 1 f f 

G(t,x, z;t' ,x , z) = —9(t — t') ——Or / dm J 



o 



(2tt) 2 Vm 2 + k 2 



xe ikR J 2 (mz)J 2 (mz'), R=\x-x\. (6.2) 
Taking into account ( p. 7| ) we find that R from the previous equation satisfies 



R 2 = \ f- ?\ + (x 3 - (x 3 )') 2 (6.3) 



31 We would like to thank Robert Myers for bringing this to our attention. 
32 It should not be confused with the IR cutoff of section ||. 
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where x 1 and x 2 are the coordinates on the transverse plane and x 3 is the collision axis. Finally, 



we perform the k integration of (6.2). Working as in [92] we obtain 



G(t, x, z\ t', x' , z') 



AnR 



6{t - t')d R 9{{t - t') - R) 9{{t -t') + R)— / dm 



x m , 



J (my/(t - t') 2 - R 2 ) J 2 {mz) J 2 {mz') 



(6.4) 



It is easy to check that ( |6.4j ) is the correct Green's function using dimensional reduction: one may 
apply 27r J °° dRR (integrate on the transverse plane) on G. The integral is trivial since it is a total 
derivative. The upper limit cancels because of the 9((t — t') — R) and the final result is precisely the 
Green's function derived earlier in [92] . Now, differentiating in ( |6.4|) with respect to R we obtain 



oo 

G(t,x, z;t' ,x , z) — — 9(t — t') 9((t — t') + R) — / dmmJ^ijn z)Ji{m z) 
An z J 



-m* (my/ it-*)* -# ) m - L l- R { + m - ^ - R) 



R 



(6.5) 



where we have omitted one vanishing term. This Green's function shows explicitly that the alter- 
ation of the metric occurs for positive times and exists only in and on the forward light cone - it 
is causal as it should be! The second term on the right-hand side of ( |6.5|) does not mix x^ with z 
directions and it is localized. This term may modify the part of the metric on the light-cone, that 
is it may modify the shock waves. However, as we will see, this term does not contribute to the SE 
tensor of the gauge theory at the order fi 2 that we are considering. Any (possible) modification of 
the shock waves occurs in the bulk and not on the boundary 33 . Actually, this is one of our main 
conclusions for this paper and we will analyze it extensively in section 0. 

Equation ( p.5|) provides the most useful form of the retarded Green's function for our pur- 
pose. Nevertheless, integration over m is possible and so we perform the integral for completeness 
obtaining 



G(t, x, z]t , x, z) 



9(t - t')d R [0((t - t') - R) 9{{t - + R) ^0(s)9(2 - *)^7==S 

(6.6) 



with 



3 The modification would not come with a z 4 dependence but with a higher power in z. 
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{t - t'f - R 



,z — z 



l\2 



2zz' 



(6.7) 



The Green's function of ( |6.5|) will be in general convoluted with powers of z. These integrals 
over z may be performed and as we will see, the results involve hypergeometric functions that 
terminate when the power of z is even. We have 



r°° i 

J z ,n G(z; z')dz' = —9(t- t') 9((t - t') + R) 

x 9((t - t') - R)n n ((t - t' f - R\z 2 ) + z 



,2 D 2 2^ , J((t-t')-R) 



R 



1 

47T 



9 (a; + — 9 [x — x' ) 9 i y2{x + — x' + )(x — x' ) — |f — r' 



x n n (2(x+ - x /+ )(a;- - a;'") - |f - rf , z 2 ) + f) ^ 



R 



n = 6,8,10 



(6.8) 



where in the second equality we have used (|6.3| ) while we have defined 



Il n (u 2 ,z 2 )=(n-4)u~ 3+n / 2 z 2 



( 2 2- 
[U - Z 



F(2-n/2,3-n/2,l,z 2 /u 2 ) 



+ (-u 2 + (n- 2)z 2 )F(2 - n/2, 3 - n/2, 2, z 2 /« 2 ) 

In equation (|6.9|) m 2 is defined by 

m 2 = 2(x + — x' + )(x~ — x'~) — \r — r'\ 2 



n = 6,8,10. (6.9) 



(6.10) 



where f is given in ( |2.7| ) while F denotes hypergeometric functions. The above equation may be 
simplified in the case where n is even and reduces to a simple polynomial which is not necessary to 
display explicitly here. Instead, we will restrict just to the z 4 coefficient since according to ( p.ll|) 
it encodes all the information (of the SE tensor of the gauge theory) we want. Expanding in z we 
obtain 



nj 4 = -n(n - 2)(n - 4) [2(x + - x' + )(x~ - x'~) - \r - r'\ 2 



-3+n/2 



(6.11) 
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Indeed, we see that for even n, (|6.11 ) is a simple polynomial of x M . The second piece of (|6.8|) shows 



> 4) dependence and so according to ( |2.11|) does not contribute (to the 



second order in /i) to the gauge theory SE tensor. However, it does modify the metric in the bulk 
and as a result it contributes to T^ v at higher orders in fi. Having performed the z integration, we 
proceed to the rest of the coordinates. We organize the remaining integrations in the following two 
subsections but we will leave the details for the appendices. 

6.2 Integration over the Light-Cone Plane 

We begin by introducing the following compact notation 

-I poo /*oo /* 

G ® z' n f(x'n\ zi = dx l+ J dx'- J dV (U n (x K - x' K , z)\ z4 ) f(x'n 

x 9 (x + - x' + ) 9 [x~ - x'~) 9 (\/2(x + - x'+)(x- - x'~) -\r- r'\ \ , (6.12) 

where /(x M ) is any arbitrary function of x^ while the last integral denotes integration in the trans- 
verse plane. We wish to perform the x integrations for all the possible cases that we will encounter 
while specifying g^ v from ( |5.8|) . We organize these integrations in five cases while we leave the 
details of the calculation for Appendix 

Remark 1: 



The results of all of the integrations in the light-cone plane (performed in appendix [A]) are 
proportional to the product 9(x + )9(x~). This implies that the second order corrections to g^ v 
appear in the forward light-cone which is what we had initially demanded by seeking a causal 
solution (see ( |2.13| )). 



Remark 2: 

The right-hand sides of ( |5.8| ) contains expressions of the form t\t<i differentiated with respect 
to x^ in some fashion. According to fl2.13| ), these expressions are proportional to 5(x + )5(x~) 
or their derivatives. Our previous analysis 34 has already taken care of the x ± integrations and 
so beginning from the following subsection, by we will mean just the transverse part of £1^: 

-n hg(^(x i ±by + (x 2 ) 2 ) . 

6.3 Integration over the Transverse Plane 

Having performed the x integrations we move to the integration over the transverse plane. The 
quantities we have to integrate have the structure 35 (^(t^) or d 2 i ^(£1^2) where a,c;i,j = 1,2. 

x a x c 

The subscript a (c) and the superscript i (j) denotes differentiation of the source t a (£ c ) with respect 



34 Along with the work of Appendix (|X|). 

35 There is another case where we have to integrate terms of the form V^ti^- However, according to (4.12) these 
(transverse) integrations arc trivial as they involve delta functions. 
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to the space-time coordinate x l (x^). We may reduce the number of the different integrals that we 
have to perform by working as follows. We firstly introduce the vectors 



61 = (&11A2) b 2 = (b 21 ,b 22 ) (6.13) 
and generalize the form of (the transverse part of) t\p given by ( |2.13 ) to 



ti(r — bi) = —fj,\og(kr 1 ) t 2 (f — b 2 ) = -fi\og(kr 2 ) (6.14) 



where fx }2 were defined by ( |4.18| ) . The next step is to exchange the derivatives acting on ti t2 , that 



is d x i with differentiations with respect to b's of ( p.13 ), that is with —db ai 36 . Finally, at the end of 
our calculations we take the limits 

£ -*•(&,()) b 2 ^{-b,0). (6.15) 



Looking at the right-hand side of equations fl5.8| ) we see that they involve the product t\t 2 differ- 
entiated with respect to the transverse coordinates 37 . Exchanging the transverse differentiations, 
according to our earlier discussion in this subsection, with derivatives with respect to the compo- 
nents of b\$, we see at once that we have to calculate the following integral 



-I /*00 /*27T 

M (n) (n,r 2 ,r) = — / / r'dr'd(j>'e{T - r') | (r 2 - r' 2 ) n log(Jfe|P + rl |)log(Jfe|P + f 2 \) (6.16) 
27r Jo Jo 

with n = 0,2 38 and where we have introduced the convenient factor 77-. The integrations are 
performed in Appendix [B] and the final result has the form 

M {n \ ri ,r 2 ,r) = 9(n - r)6(r 2 - r)A^ n) (r u r 2 ,r) + 9(r - r 2 )9{n - t)M { 2 ] \r u r 2 ,T) 

+ 0( T - ri )6(r 2 - r)Mf\ ri ,r 2 , r) + 6(r - r x )e{r - n)M^\n,r 2 , r) (6.17) 

where the Aif^'s may be found using flC7|) , (|C8|) , (|C^ ) and (|C10| ). This is the last ingredient we 
need which allows us to obtain the desired solutions for equations ( J5.8|) . We display the results in 
the next section. 



36 So for instance t l x it 2yX i takes the form d% b (tib^)- 

37 Where £1*2 ~ log(fcri) log(fer2) (see ( |2.13 ) while the a; ± contributions have already been taken into account in 
the previous subsection (see Appendix |A|) . 

38 The superscript ( n ' on My 1 ' does not denote the order in the expansion in powers of /1. 



- 29 - 



The final formula for T, 



Calculating the SE tensor is the main goal of our paper. The basic formula for our purpose is 
equation (|2.9| ) whose right-hand side is specified using ( j5.8|) . The necessary equations we need are 
equations O, fH), Q), (HI, (HD, ®, ®, O, S), (PD> (ED, along with 
( |C8|) and ( p9| ) along with QC10 ). Since we are interested in the z A coefficient of only, our starting 
point for specifying any component of is equation (|6.12|) . In particular for T\\ we use ( |A2a| ) 
along with (|C7|) and and also (|A"2bD along with Q and QCTCj) . For T 12 we need along 
with ( |C9| ) and ( pidj) . For the T ++ component we need (|A6|) along with (|C3|), ( |A4b| ) along with 
an d also flSgp . Finally, for T +1 we need QSJJ) along with (£7|) and Q and ( [MaD along with 
( [4.12| ). At the end of the integrations we should take the limits as in ( |6.15| ). The final result is 



n: 



2tt 



6i, 2 ->(±6,0) 



TV 2 

2 iV c / -\2/ 



2tt 



2J(r) + 3 



fr 2 - r 2 ) 2 



r 1 



r - n : 



X 



&i, 2 -K±b,0) 

(fen - fegiK^ - &n) + (612 ~ MO^ 2 - &12) 
|fci-& 2 | 2 



(Tj'V)) = 2 y u 2 ^%x-^(s + )0(a;-) lim <J 4 (<9 b21 - 2d 6u ) M (0) + 3^ 



n: 



r? — r 2 



2tt 2 



6i, 2 -)-(±6,0) 



r - n 



x 



(&11 - &2l)^ 



r? — r 2 



&i-fe 2 | 2 



(6.18a) 



<Tf(^)>=4/^(*+)0(^ lim [2^ 22 +2C fe21 -C. 12 -CJ-^ (2) > (6-18b) 



(6.18c) 



+ 2(6n - a: 1 ) log (fcl&i - 6 2 |) J, (6.18d) 



(6.18e) 



where J = J(t, r 1; r 2 ), -M (0) = M^(r, n, r 2 ) and M {2) = .M (2) (t, ri, r 2 ) are given by (|C3|), 
8h and (|U9|) respectively while the equation for T+_ follows from the tracelessness condition, 



equation ( |5.3| ). The components T^ 2 are obtained from TV when interchanging the labels 1 f>2 

(2) (2) 7* ~* 

while components are obtained from when interchanging simultaneously the +,61 -H- —,62 



before taking the limits &i j2 — > (±6,0). Formula ( |6.18| ) is our final result which we will analyze in 
the next section. 
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7. Discussion 



7.1 Regime of Validity 

In section |3| we argued that the SE tensor of the single nucleus can be trusted for transverse distances 
r < \. Therefore, in the case where we have two nuclei colliding with an impact parameter b, we 
trust the solution for 



r, r li2 «i r li2 = T b) + (x 2 ) 2 . (7.1) 

On the other hand, from dimensional analysis, we expect that higher powers in \i should be 
compensated by higher powers in r, t\ and r 2 which generally will be multiplied by logarithms with 
arguments &?"i,2, kb and kr which we will collectively call q. In particular, we expect that 

the SE tensor at mid-rapidity should behave like 



T,u = /i 2 £ c% mnp T k b l r?r" 2 \ogP[q<§ mnp ] + // £ 41n P "^« ktffeSLj + . . . (7.2) 

klmnp klmnp 

where k, I, m, n are integers and Ckimnp real numbers. In each sum the superscript (j) where 
j = 2, 3, . . . denotes the order in the expansion in powers of /i. From dimensional analysis, at any 
given order in /i 39 , the powers should add up to 3j — 4, i.e. k + I + m + n = 3j — 4. In addition, 
in order to have a finite energy production at any finite piece of the space-time, we expect that 
when a negative power of one of the variables r, r± or r 2 appears, it should be compensated by 
a different variable with positive power that tends to zero faster; the term in the square bracket 
of equation (|6.18c|) provides such an example. Mathematically this is achieved by the presence of 
theta functions. As a result, this fact shows how causality and conservation conspire in order to 
prevent the produced T uu from becoming infinite. However, this does not mean that r\ and r 2 
cannot be larger that r. What it is implied here is that when r 12 > r then r 12 will appear as 
denominators while r as a numerator; the term (8(r 2 — T\)Q(r\ — r) + Q(t\ — r 2 )6(r 2 — r)) M.f\r) 
of equation with (r) given by ( |C7a| ) is an example of this case. 

In analogy to [93] and [92], we may imagine the gravitational description of the process as 
a shower of graviton exchanges between the sources and the bulk where the gravitational field is 
measured. Below, we will argue that an expansion of the form of equation (772) is consistent with 
our diagrammatical intuition. The first term of the right-hand side of equation ( [7.2[ ) corresponds 
to figure ||. At central rapidities the effective vertices (see figure [l]), loosely speaking, should have 
as dimensionless couplings the quantities 



HT 3 \log(k ri )\ /ir 3 |log(fcr 2 )|. (7.3) 



:!'.) 



Recall that fi has dimension of mass cubed. 
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We firstly note that the couplings of ( \l.3\ ) being small is equivalent to 



fir 3 \ log(Ari)| < 1 fir 3 \ log(Ar 2 )| < 1 /ur 3 | log(fcr)| < 1 fir 3 \ log(kb)\ < 1. (7.4) 

We stress that the last two terms of the previous equation are a consequence of the first two and 
are not extra conditions. Now, because of ( |7.1| ), the absolute value of the logarithms is greater than 
unity and result (\IA\) holds provided that 



/2T S < 1. (7.5) 

Equation ( [7.5| ) justifies most of the terms in the expansion of ( |7.2| ). The negative powers of r 12 and 
r exist by dimensional analysis while their finite contribution to T^ v has been already justified by 
causality and conservation. 

We still have to investigate the conditions under which the inverse powers of b (see (|6.18c| ) 



and ( |6.18d|) ) do not cause the breakdown of the perturbative expansion (in ji). It seems that the 
inverse powers of b appear when the logarithms are differentiated with respect to x 1 or x 2 and are 
compensated by x 1,2 6'(r — r 12 ) or r — r 12 #(r — r 1>2 ). Thus that as long as 



/ir 3 ^ < 1 (7.6) 



we may trust the first terms of (|7.2| ). On the other hand, unless we calculate the next order 
correction of T^ u in /i, we cannot be certain whether the inverse power of b appearing in (|6.18c| ) 



and (|6.18d|) iterates or whether is a harmless overall factor appearing just once at order fi 2 . In the 



latter case, restriction (|7.6j ) is not required. 

As b — > the metric exhibits in addition a logarithmic divergence. From a physics point of 
view, a similar logarithmic divergence is observed in gluon exchanges between two sources (such as 
quarks) in the eikonal approximation as the two sources approach each other [174]. Both types of 
divergences are due to the corrections of the bulk SE, Jmn- 

It is also interesting to consider the high energy limit. Taking into account ( |3.17| ) and ( ^.19[ ) 
we have that 



^ > k. (7.7) 
We note that the inequalities f[7.5|) and (|7.1|) are consistent with (|7.7| ) and when combined yield to 



fcr</isr<l. (7.1 



-32- 



What one has to keep from the above discussion is that the first terms of T^ u in ( |7.2|) hold for 
sufficiently small proper times and high energies (compared to the transverse scale k) such that 
inequalities (|7.1|), ([7.4|), ( |7.6| ) and ( |7.7| ) apply simultaneously (equivalently when ( |7.8| ) applies). In 
particular, (7A_) justifies the choice of the transverse profile of the and ( |7.7| ) the high energy limit. 



Finally, ( |7.4| ) associates the linear approximation of the gravitational field with a small proper time 
expansion. This makes physical sense since the gravitational field in the linear approximation is 
initially good enough to describe the process until the nonlinearities of the field will dominate for 
later times; this is just like the McLerran-Venogopolan model [29] when analyzed at early times [147]. 

7.2 Investigating the solution 

The energy density e of the produced medium may be found by the formula 

e = Too = \ (T++ + 2T + _ + T__) = i (T++ + T n + T 22 + T__) (7.9) 

where the second equality is the tracelessness condition (|5.3|). Unfortunately a complete analysis of 
the results appears to be rather involved and is left for a separate project. Nevertheless there are 
some kinematical regions in which we can derive some conclusions. In particular, we investigate the 
energy density e for the kinemtical regions r «C rx,r% and ri,r 2 <C r. According to figure ^ these 
two regions belong to the regions / (/') and III (III ) respectively 40 . Certainly in either case, 
Tj), ([7.4|) , ( |7.6|) and ( [7.7|) are assumed to be still applicable. 



We argue that the leading contributions to T^ v in the two regions that we will investigate, come 
from the quantities J and A^ -* which arise as a consequence of the initial shock waves and 
not because of Jmn or A4^. The reason is because these quantities come with a double (large) 
logarithmic enhancement in the kinematical areas under consideration. The rest of the terms in 
T^u either vanish exactly because of the presence of 9 functions or do not exhibit this enhancement. 
Table ^exhibits the double logarithmic behavior of Ai^. An analogous enhancement exists for the 
factor J given by equation (|C3a| ) . The logarithmic suppression of the quantity Ai^ 2 \ equation (|C9|) , 



is not evident. Thus, we have constructed table |3| and table |], in appendix || in order to clarify the 
behavior of Ai^ 2 \ One has to keep in mind that when appears, it has to be differentiated twice 
as shown in ( |6.18a| ) and ( S.18b|) and then appropriate limits are taken as in ( |6.15| ). This procedure 



lowers the logarithmic enhancement of these terms so in this approximation they are ignored. 

7.2.1 Energy density for regions / and / 

We are interested in the energy density at the regions defined by 41 



(7.10) 



40 The reason they come in pairs is because there is an obvious r\ ■<-> ri mirror symmetry. 
41 Equation (7.1) still applies. 
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Figure 6: The reaction plane: Region III corresponds to 7*1 < r 2 < r while region III' corresponds 
to r 2 < n < r. The dark dots are the centers of the nuclei at impact parameter 2b while ri >2 denote 
the distance of the arbitrary point r from the center of each nucleus (right and left respectively). At any 
given proper time r, the propagation from the centers will reach the points on the peripheries (at most). 
This suggests that any given point f on the transverse plane of the produced medium at given r will 
evolve according to the region to which it belongs. We find that our calculations are consistent with this 
expectation (see (|7.25|) ). 

According to figure [| we are looking at (a sub-region of) regions I and I . Mathematically this 
means that we are interested in the terms that are multiplied by 9(rip — r). 
T11 and T 22 contributions: 

The largest contribution to Tn (and T22), according to our previous discussion and equation 
( pM8a| ), comes from the coefficient Ai^ . Using table |2] we find 

N 2 r 

Tn = T22 = (0(r 2 - n)e(n - r) + e(n - r 2 )#(r 2 - r)) r 2 log(Axi) log(A;r 2 ) 

+ 0(/i 2 ^)}. (7.11) 
T ++ and T contributions: 

In this case the term that gives the largest contribution to T ++ is the term proportional to J 
of appendix 0. This term contains the factor log(A;ri) log(A;r2) and hence we deduce that 
42 Similarly for T__ whose main contribution is obtained from that of T++ by —6, -h- b and + -H- — . 
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N 2 r 

T±± = -2// 2 -f-(x±) 2 | (0(r 2 - n)0(n - r) + 0(n - r 2 )6{r 2 - r)) log(fcn) log(Ax 2 ) 



r x r 2 ' 



(7.12) 



The last step is to substitute ( |7.11| ) and ( |7.12| ) in ( |7.9[) to find that 



7V 2 r 
e 7 = = V^f(2r 2 - {x + f - (x_) 2 ){ (0(r 2 - n)0{n ~ r) + 0(n - r 2 )#(r 2 - r)) log(fcn) log(Ax 2 ) 



+ 0(/* 2 



rir 2 



(7.13) 



We observe that the energy density is symmetric under 7*1 -H- r 2 in regions 7 and 7 . The diagonal 
elements of T^ v in the 3 coordinates at central rapidities where (x + ) 2 ~ (x ) 2 

obey 



±r 2 



00 



22 
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33 



(7.14) 



and receive their main contribution from the center of each nucleus separately in a multiplicative 
way as equations (|7.11 ) and ( [7.12 ) show. 



7.2.2 Energy density for regions 777 and 777 

Analogously to the previous case we are interested in the energy density at the regions defined by 43 



T\ t r 2 T 



(7.15) 



which implies 



x , x 2 , b <^ T. 



(7.16) 



Inequality ( [7.15 ) implies that the propagation from the center of both nuclei has enough proper 
time r in order to travel and affect the arbitrary point r on the transverse plane. The inequality 
( |7.16|) shows that for the case under consideration, we may derive results for the area around the 
origin, that is around the center of the collision. This kinematical area concerns regions 777 and 



Equation (7.1) still applies. 
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Ill' of figure |6]. Mathematically this means that the contributing terms should be multiplied by 
0(r-r 1)3 ). 

Tn and T 22 contributions: 

Working as we did for regions / and I , we find that 
N 2 r 

T n = T 22 = ^ 2 -^{ (0(t - r 2 )6(r 2 - n) + 0(r - n)0(n - r 2 )) r 2 log(Arr) log(fcr) 

+ 0(//Vlog(Ar 1 , 2 ))) (7.17) 



where we have assumed that | log(A;ri j2 )| ~ | \og(kb)\ <C | log 2 (kr)\ although in this region | log(A;ri i2 )| 
| \og{kb)\ > | log(fcr)|. 

T ++ and T contributions: 

The procedure is the same as before. Employing equation QC3a|) we have 
N 2 ( 

T±± = -fi 2 ^(x±) 2 [ (6(t - r 2 )6(r 2 - n) + 6{r - n)e{n - r 2 )) log(A;r) \og{kr) 

+ 0^ 2 —)}. (7.18) 
nr 2 ) 



Finally, substituting ( |7.18| ) and ( |7.17| ) in equation ( [7.91 ) we deduce that 



N 2 

em = e IIV = fi 2 ^(2r 2 - (x + ) 2 - (x_) 2 )\ (B(r - r 2 )6{r 2 - r a ) + 6{r - n)0(n - r 2 )) log(fcr) log(fcr) 



7T 2 



+ 0( /U 2 r 2 log(fcr 1 , 2 ))}. (7.19) 

The conclusions are similar as in the previous case. The difference is that the arguments of the 
logarithms (nuclear profiles) have been changed to proper time r. This means that there is enough 
proper time r for the propagation from the center of each nucleus to the observation point. At 
central rapidities we again have that 



Too = T n = T 22 = -3T 33 . (7.20) 
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7.3 Momentum anisotropy and spatial eccentricity 

The momentum anisotropy e p (r) and spatial eccentricity e x (r) of the produced medium are defined 
by 



rdxHx\T n -T, 



221 



f™ oo dx 1 dx 2 (T u + T 22 ) 



IT 



dx 1 dx 2 T 00 {{x 



1\2 



/^dxWTooCOr 1 ) 2 + (* 2 ) 2 )' 



(7.21) 



Having found the formula of T^ u , we are in the position to plot both of these quantities as a function 
of proper time r and the impact parameter b. As the underlying expressions are complicated, we 
leave this for a future paper. For the moment we restrict ourselves to finding the momentum 
anisotropy and spatial eccentricity in the two extreme cases discussed in the previous subsection. 

The equality T 1X = T 22 in both of these cases applies very close to the collision center or very far 
from it. Both are considered at sufficiently early times (compared to the energy of the collision/see 
subsection |7. 1| ) . Formula ( |7.21[ ) shows that that the momentum anisotropy in either case is zero, 
that is 



(t) 







6>r 



(r) 



0. 



(7.22) 



The physical meaning of this result is that at the early stages of the collision (where our analysis is 
applicable) and at the center of the collision, the system behaves isotropically until it will starts to 
feel the anisotropy from the edges of the nuclei at later times. This result agrees with [1] obtained 
using hydrodynamical methods. On the other hand, applying the approximate formula ( |7.19| ) for 
the energy density in equation (|7.21| ), we may estimate the eccentricity e x for regions /// and III' 
of subsection |772|. We find 



IL 



dx'dx^ix 1 ) 2 - (x 2 ) 2 ) 



b<7 



IL 



DmUDjjj 



i dx l dx 2 ((x 1 ) 2 + (x 2 ) 2 ) 



(7.23) 



where Dm U Djjjr is the two dimensional surface defined in figure |] as the intersection of the circles 
of radius r (almond shape). Thus the integration area is r dependent. In principle, we should 
add the contributions to e x from the other regions. Unfortunately, performing the corresponding 
integrations (using T 00 from ( |6.18| )) analytically is hard and so is left for a future project. We argue 
that in a realistic situation of a heavy ion collision the important contribution to spatial eccentricity 
comes from the overlapping region (the almond) while whatever is left out is less important. Hence, 
equation (|7.3| ) may suffice as a first approximation in determining e x . Assuming this and performing 
the integrations of (|7.3|) and employing formula ( [7.1 9|) we find that 



2 Vx 2 - 1(1 + 2s 2 ) - 3a; 5 



sec 



[x 



6<T 



3x 2 



-3VX 2 



1 + (2 + x 2 ) sec" 



x 



T 



(7.24) 
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Figure 7: The left plot depicts the eccentricity e x given by equation (7.23) as a function of x = ? for 
x ^> 1 (we assumed r > 106). It approaches asymptotically to zero at large x. The right plot depicts e x 
(also from equation (7.23)) as a function of ^ = | for (again) i>1 (we assumed ^ < 0.5). It provides 
the same information as the plot on the left and it implies that for fixed r and as b — > 0, we have that 
e x —> 0. In both figures, the dashed lines represent the asymptotic behavior of e x at regions where our 
approximations are no longer true as either r is not large enough (compared to b and hence (7.23) is not 
applicable) or because inequality (7.6) is violated (large r). 



Equation ( [7.24Q shows that e x is a function of the ratio x = | while the result ( [7.24Q applies for 
x ^> 1. In figure ^ we plot e x as expressed in ( |7.24|) as a function of x and also as a function of - for 
x > 1. For smaller values of x equation (|7.24 ) is not applicable. Nonetheless, this constraint can 
be easily removed and our result may be extended to smaller values of x (that is r) by using the 
full result for T^ u given by equation ( |6.18| ). This is left for future research (due to the complexity 



of ( p. 18 )). From the figure we see that as r increases significantly compared to b (left figure) or 
equivalently as b — > for fixed r (right figure), we have that e x — > 0. The physics behind this is 
that when r dominates over b, the eccentricity e x is insensitive to b which expresses the asymmetry 
of the system and hence e x tends to zero. On the other hand as x moves towards even larger values 
(r increases), equation ( |7.23|) is no longer valid because linear gravity is not applicable any more 
(that is inequality ( |7.5|) is violated). In order to overcome this constraint, one needs to include 
higher order terms in /i in our perturbative solution. Nevertheless, figure |7| describes qualitatively 
eccentricity for these (intermediate) values of x = t and it shares the same asymptotic behavior 
with results already obtained in the literature using pQCD methods [150, 151] and hydrodynamic 
simulations [1], [152,153]. 



7.4 Conclusions 

Using the AdS/CFT correspondence we predict the evolution of the SE tensor of the produced 
medium in heavy ion collisions taking into account a nontrivial transverse coordinate dependence. 
Our predictions apply for sufficiently early (see subsection |7. 1|) proper times. Our conclusions are 
summarized as follows. 
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1. The SE tensor T^ v (to 0(/i 2 )) is given by equation ( |6.18| ) and evolves nontrivially constrained 
by causality in an intuitive way. In particular, the behavior of T M „ at any point on the transverse 
plane, is determined whether the propagation from the center of each individual nucleus has enough 
proper time to reach the point under consideration or not. Figure || is a snapshot taken at given 
proper time r and depicts the six kinematical regions where T MJ , evolves differently while it has the 
general form 

TjS = fi 2 6(x + )6(x~){6(r 2 - r a )0(n - t)A^{x*, b) + 6(r 2 - t)6(t - n)Aji(x K , b) 

+ 6(r 2 -r 1 )6(T-r 2 )A I JJ(x K ,b) + (b^ -b)}. (7.25) 

The indices J, II, III on A^ u correspond to the regions J, II, III of figure || respectively 44 . 

2. Our early (proper) time analysis (see subsection ( |7.1| )) is not adequate to prove or disprove 
wether a thermal medium is produced allowing for its hydrodynamical description. In principle, 
one could get a hint for the formation of a horizon and hence for thermalization by plotting the 
(components of the) metric cjmn up to the order /i 2 , which is what has been computed in this work 

45 , and investigate the possibility for the metric to change its sign on some (generally) hypersurface 

46 . If this is the case, then the set of points on this hypersurface would correspond to some 
temperature of the medium which would generally depend on and b. Such a procedure could 
give an estimation for the thermalization time; this is non-trivial because of the complication of the 
corresponding formulas and it is therefore left for future investigation. However, even if a change 
in the sign of the metric appears, this wouldn't suffice to prove thermalization but it would merely 
be a good indication that it occurs. This is because of the incomplete analysis, that is the linear 
approximation, which we have used in the present work. In order to have a complete proof one 
should solve the problem to all orders in /i, that is the full non-linear problem. 

3. The energy density simplifies at r t ^> r, i—1,2 and r, <C r where is the distance of the 
arbitrary point r*on the transverse plane, from the center of the nucleus i where i=l,2. It behaves 
in a way proportional to the product of the SE tensor of the two nuclei. In particular for these 
two regions (r^ ^> r and <C r) we get that the energy density evolves as in ( |7.13| ) and ( [7.19| ) 
respectively. More specifically, the energy density for r has a logarithmic dependence with r. 
Our result has the same log 2 (A;r) dependence as [147] (see also [148,149]) which has been derived 

44 The (6 -f> —b) terms cover the rest three remaining I', IV and III' since under this interchange ri f-> ■ 

45 Actually in this paper we have only computed the z 4 coefficient of guN to 0(/j, 2 ). Computing the rest powers 

in z (to 0(fi 2 )) may be achieved similarly. 

46 This hypersurface, if exists, may be in principle found by an analytic computation. Plotting the metric would 

be a convenient way to locate the hypersurface but it is not a necessary step. 
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using pQCD. However, our result (at central rapidities) has in addition an overall r 2 dependence. 
Thus one may conclude that (in the aforementioned kinematical window and in the limit b — > 0) 
the energy density, except from an overall r 2 dependence, is symmetric under r -H- r. In fact, this 
conclusion is applicable to all the regions where our analytic formula for T^, equation ( [7.25Q , is 
applicable. A similar behavior has been observed in [108]. 



4. The formula we found for T^ v (see (|6.18|) ) and which is valid at early proper times allows us 
to compute the momentum anisotropy and spatial eccentricity as a function of the ratio In this 
paper, we perform the calculations for both quantities for r ^> b 47 . We find that our results (see 
( |7.22|) , ( |7.24j ) and figure |7|) agree qualitatively with results in the literature obtained in the context 
of pQCD [150,151] and hydrodynamic simulations [1], [152,153]. 



5. We have seen that the SE tensor of the produced medium depends strongly on the product 
of the SE tensor of each nucleus before the collision. What is more interesting is that when the 
proper time is large compared to the observational point, the argument of the profile is replaced 
by r. We conjecture that this conclusion is independent of the choice of the transverse profiles 
under the logical hypothesis that the nuclear density is greater in the center of the nucleus than the 
edges. Assuming this in a more realistic situation, one may take the Woods-Saxon energy density 
of a nucleus at rest given by 



Too = e ws (r, x 3 ) « M mE , x = (x\x\ x 3 ) (7.26) 

1 + e a 

where a and R are positive constants and M has mass dimension four (compare with fi that we used 
in this paper) while the remaining components of T^ u are zero. The next step is to boost this tensor 
(see also [114], [140] and [142]) into the infinite momentum frame keeping the ratio \i 



where , 1 is the boost factor, fixed. Assuming the boost is along a; 3 , the result of such a boo: 
gives 



OG 



T_ = /!*(*-) / dx 3 VrW)2 - , r 2 = (x 1 ) 2 + (x 2 ) 2 . (7.27) 

J -°° 1 + e s 

Consequently, according to our conjecture, we believe that at early times the energy density 
should depend strongly on the factor 

/OO POO 
dx 3 e ws (i] > ,x 3 ) / dx 3 e ws (^>,x 3 ), £> = max{r 1 ,r}, 7?> = max{r 2 ,r} (7.28) 
oo J — oo 



47 



But still sufficiently at small r. 
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where r» is the distance of the observational point from the center of the nucleus i (i=l,2). Hence at 
the center of the collision where the proper time r satisfies r > r 12 , we conjecture that the energy 
density (and generally T^ u ) behaves like 



e t 



flT 



dx" 



1 + e 



1 

y/T 2 + (x 3 )' 2 -R 



(7.29) 



6. Any terms in equation ( |6.18|) that give not proportional to J , Ai^ and Ai^ are due to 
Jmn- As is evident from these terms, the presence of the bulk sources Jmn and the back- reactions 
affect g^v and hence not only on the forward light-cone but also inside. This implies that we 
cannot in principle solve Einstein's equations in vacuum arguing that we are away from the sources 
unless we know the boundary conditions that these sources enforce on the metric. This is in analogy 
to classical electrodynamics: solving Laplace equation for the scalar potential away from a point 
charge sitting at the origin without specifying the boundary conditions, one may obtain the trivial 
(zero) solution which obviously is not the correct one. 



7. The presence of the impact parameter b is a necessary requirement and not an additional 
complication introduced in the problem. This is obvious from the solution of the SE, equation 
Q3.181 ), tensor which diverges in the limit b — > 0. More importantly in the same limit, according 
to (|4.20|) , the conservation of Jmn is violated. This suggests that a head-on collision may not 
be investigated using classical gravity. Instead, one has to apply a quantum theory of gravity in 
the same way one can not predict the electron-positron annihilation (in head-on collisions) using 
Maxwell's equations. One has to turn to Quantum Electrodynamics in order to describe the process. 



8. We still face the problem of negative energy densities appearing in [90] - [94]. This is implied 
by the fact that while conservation (see ( |5.7|) ) applies, the initial shock waves continue to propagate 
in the forward light cone unaltered. Since these carry the total energy of the collision, we conclude 
that the net energy production should add to zero. This should not be a surprise but rather a feature 
of our approximation; it is just a consequence of the eikonal approximation which we consider in 
this work. Since the eikonal approximation is always valid at early times (see [174] for some eikonal 
results), negative energy is always present. 



9. For future projects we propose that one could compute the spatial eccentricity and momen- 
tum anisotropy at early times extending our preliminary results (see conclusion 3) by making use 
of the (full) analytic formula ( |6.18| ) that computes T^ v (for sufficiently early times). This would 
extend the result of the left plot of figure |7] to smaller values of x = | (that is to the left). In 
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addition, a similar analysis could be applied to the shock waves that correspond to (|3.9|) (see [114] 
and [115]). These shock waves have the advantage that they are created by a point-like SE bulk 
tensor and hence the back-reaction should be specified by the geodesies. Unfortunately, their form 
is very complicated and the analysis would be much more involved. 

10. Finally, our technique has already been applied to shock wave collisions in ordinary four 
dimensional gravity [175] , taking into account back- reactions. The problem of axisymmetric colli- 
sions in four dimensional gravity has been analyzed extensively while the solutions were obtained in 
vacuum [137-139]. In [175] one may see the way solution of [137-139] gets modified in the presence 
of matter and non-zero impact parameter. 
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A. Solving equations (5.8) 

In this appendix we perform the x ± part of the integrations resulting when the Green's function 
acts on the right-hand side of (|5.8| ). Defining 



p = \ f' — f\ (Al) 
and proceeding as in (|6.12|) we find out that we have to deal with five different cases: 



Case I: z n 6(x + )S(x ) terms 
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G ® z' 6 5(x' + )5(x'-)f(f') = 3z 4 9(x + )9(x-)— [ d 2 r9{r - p)f(r'), (A2a) 

2ir J 

G <g) z'*5{x' + )5{x'-)f(r') = 12z 4 9(x + )9(x-)^- [ d 2 f9(r - p)(r 2 - p 2 )f{f). (A2b) 

27T J 



Case II: z 8 5'(x + )<5(x ) terms 



G <g> z'^V+WO/^) = 24aT,2 4 #(:r + )#(ar)-^- / d 2 f<9(r - p)/(f"), (A3) 

2n J 

where we have integrated by parts the S(x' + ) and have dropped two vanishing terms. 
Case HI: z n (x-) m 5(x + )9(x-) terms 

G ® z' 8 5(x' + )9(x'-)f(f') = -%z A 9(x + )9(x-) [ d 2 r9(r - p)(r 2 - p 2 ) 2 f(f') (A4a) 

n t z J 

G <g> z m x'-9(x'-)5(x + )f(r') = l^Q-z 4 9(x + )9(x-)^- [ d 2 f9(r - p)(r 2 - p 2 ) 2 f(f') (A4b) 

2 r 4 2tt J 

where in both cases we have used the identity 

/oo 
dx'~9(x- - x?-)6(y/2x+(x- -x'-) - p)9{x'-)f(x'-) 
-oo 

= 9(x + )9(x~) [ dx , -9( v / 2x+(x- -x'-)- p)f(x'-) 
Jo 

= 9{x + )9{x-)9{t - p) / dx'-f{x'-) (A5) 

Jo 

Case IV: z 8 S"(x + )S(x~) terms 



G ® z' s 5"(x' + )5(x'-)f(f r ) = ^z 4 9(x-) (d 2 x+ ) I J°° dx' + J d 2 f'5(x' + )9 (x + - x l+ ) 

x 9( y / 2(x+ - x'+)x- - p) (2(x + - x' + )x~ - p 2 ) J/(f") 

= ^z 4 9(x-)(2x-)(d x+ ) {9(x + ) J d 2 r9(V2x~+x~ - p 2 )}/(f') 

= 24(x-) 2 z 4 9(x + )9(x~)J- [ d 2 r'5(r - p)f(f r ). (A6) 

2vrr J 
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A few explanations are in order: In the first equality we integrated by parts twice the 5-function 
and exchanged d x >+ <-> — d x + (these act on the curly bracket), in the second equality we performed 
the x' + integration and the first derivative with respect to x + ignoring two vanishing terms and in 
the third equality we performed the second differentiation ignoring one more vanishing term. 



Case V: z 8 x-6(x-)5(x + ) (t 2>x i V 2 J 1)X i + t 2 ^\t x , x 



As mentioned earlier in subsection |6.3| the expressions t\ and t 2 carry only the transverse 
dependence. We have 



IX T 4 



G <g) z' 8 x' 6{x )5(x' + ) {t 2x n V±ti iX 'i + t2 yX i2V 2 i ti yX i2^ 

= ^-z 4 6(x + ) J dx'~ J d 2 r'x'-Q{x'-)B{x- - x'~)B{ v / 2x+(x- - x'~) - p) 

X (2X + (X- - X'-) - p 2 ) (t 2<x nV 2 ± t hx n + t^aV^t^) 

'-6{x + )6{x~) J d 2 r'9 (t - ^ {x 1 - x' 1 ) 2 + {x 2 - x' 2 ) 2 ^j 

x ( r 2 _ ^ _ ^1)2 + {x 2 _ ^yy* (t 2tX nV 2 ± t hx n + t 2 ^V\t X ^) , (A7) 

where in the first equality we performed the x + integration while in the second equality we performed 
the x~ integration with the help of ( |A5|) and also we substituted ( |A1|) . In order to perform the 
transverse integrations we integrate (once) by parts the term V^ti^/i = c^/i (V^ti) with respect 
to x' 1 . I gnoring the surface term because of the ^-function, we obtain three terms. One term 
contains a ^-function term which vanishes (it has the form xS(x)). The second term is proportional 
to t 2x 'i x n Vj_£i which together with the integration by parts of the term V 2 ± t ltX i2 with respect to x' 2 
gives the contribution ~ V^V^ti. But this term, for nonzero impact parameter, is according to 
( [4.21]) zero. Hence, we are left with the term where d x a acts on (r 2 — ((x 1 — x' 1 ) 2 + (x 2 — x' 2 ) 2 )) 
and with a similar term from the action of d x a. Thus, the integrations by parts simplifies to an 
integral proportional to Vj_ti which according to ( |4.12| ) results to ~ b{f—b\) and so the integrations 
become trivial. The final result is then 

G <S> z m x'~6(x'~)S(x' + ) (t 2 yiV 2 ± t hx a + t 2 , x >iV 2 L tx, x n) 

= -12p 2 z%x + )d(x-) {x 4 )2( 1 2 ~ r } ] " ((bu - b 2l )( Xl - b n ) + (612 - b 22 )(x 2 - b X2 )) (A8) 



T 



6l 



where the fey's are defined in ( |6.13| ). 



B. Evaluating the integral (6.16) 



In this appendix we calculate the two integrals Ai^ of ( |6.17| ) for n 
the JA^ integral first 



and n — 2. We begin with 
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M {0) = / dV(0(r - \f- f'\) \og(kr[) log(£r 2 ). (CI) 

271 I 



The quantities 7*1,2 are given by Q4.18Q . The area of integration is a circle of radius r centered at 
r. The trick here is to expand the logarithms in their Fourier space: log(kr) = — J ^r^r with k 
serving as an ultraviolet cutoff. Shifting the variable r' — > r' + f, expanding both logarithms and 
performing the angular integration one obtains 48 



/ T f f rl 2 nrl 2 l piQri+ilfi | ft 

\ J \tf q*P Mr ' W+ drVj ^- (° 2 ) 

Next we perform the I and q integrations. The integrals that define J have been calculated 
in [174]. In order to perform these integrations one has to expand Jo{r'\q + l\) in an infinite sum 
of products of the form J n (r'\q\)J n (r'\l\) with n an integer and do the angular integrals (of q and 
1) first. This factors out the radial integrations over I and q into two independent integrals. Then 
one has to perform these integrations and finally sum over n. The final result reads 



j{r') = H£>k) Hv>k) + 

£>(<) 



A<V< 
£>V> 



+ Li 2 



€>V>J_ 



max(min)(ri,r') V>{<) = Tnax(min)(r2, r') 
(fi).(f* 2 ) = cos(a)rir 2 . 



(C3a) 

(C3b) 
(C3c) 



Here a is the angle between r*i and r 2 and Li 2 is the dilogarithm function. We also note that J is 
real as it should. 

We still have to integrate (|C3a| ) over r' in order to obtain Ai^. This is not trivial since 
according to the definition of £>(<) and //>(<), we have to break the r' integration into cases. In 
practice, we have to introduce 6 functions. We organize the several cases by first introducing table 
D 49 . With the help of the table and taking into account all the six possible ways one can order n, 
r 2 and r it is found that 



48 Recall from ( |4.18| ) that r- = f T b, i = 1, 2. 

49 J73 for instance means l 7(C>?7> = r'r2,^<rj < = r'rx) with J given from (|C3 



-45- 



Table 1: We defined fi,2 = (a; 1 =F b,x 2 ) 



cases 


£> 


v> 




€<V< 


Ji 


1 


r 1 


r 2 


rir 2 




Ji 


2 


n 


r' 


rVi 


r'r 2 


J2 


3 


r' 


r 2 


r'r 2 


rVi 


Js 


4 
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+ (0(r 2 - r x )Q{r x - r) + 0(n - r 2 )fl(r 2 - r)) / dr'r'Ji. 



(C4) 



The integrals involving J^, J 2 and J7a are easy to calculate. The J 4 integral is harder. It involves 
integrations of the form J d'rr'Li 2 (e ±iar ^r). In order to do such integrals we use the integral 
representation of Li 2 



Lio(z) = z 



w 



-dw. 



(C5) 



, &"-z 

Then we perform the integrals by first integrating over r' and then over the parameter w as follows 
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(C6) 



Verifying that the derivative of the right-hand side of (|C6|) is equal to the integrand is straightfor- 
ward. Now using this result, we may obtain the results of all integrations whose values in their 
indefinite form (up to a constant) are 
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( 2) 

With these integrals at hand and taking into account that M.\ (0) = 0, we write a final expression 
for M.^ . Using flC^]) we eventually have 



= #(r 2 - n)9(r - r 2 ) [^f(n) + MV2) - MVi) + Mf(r) 
+ 0(n - r 2 )0(r - n) [Ali ^) + -M^ 0) (n) - M 0) (r 2 ) + AfJV) 
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+ d (r - r 2 )e{ n - r) [M { ?{r 2 ) + M { 2 °\r) - M { 2 °\r 2 ) 
+ (0(r 2 - ri)0(n - r) + 0(n - r 2 )0(r 2 - r)) AlfV)- 



--Mf(r 2 ) 
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(C8) 



with Alf^'s given by (|C7|). 

We now proceed in evaluating the integral (6.17) for n — 2. The A^ 2 - 1 integral is calculated 
as before: one has to again break the integration into cases and so on, repeating the steps from 
equations (|Cl|) -(C8). The difference now is that the factor r' in the Al (0) integral of (|Cl|) (and 
so for the integrals of (|U7|) ) is replaced by r'(r 2 — r /2 ). Therefore, working exactly as before, one 
obtains 
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(C9) 



The .M^' s are analogous to those of flU7|) and are associated with an analogous table as table |T[ 
Their explicit forms are given by 
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In addition to the explicit form of the integrals Ai^ (equation (|C8|) ) and Ai^ (equation (|C9|)) 
that we have calculated in this Appendix, we find it useful to introduce tables ||, |] and £|. These 
tables display the behavior of Ai^ and .M^ in the six possible kinematical regions associated with 
the collision (see figure |6]) and will be useful in the analysis of the results of section [F| 
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Table 2: 
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Table 3: 
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Table 4: 



region 


ordering 


Non-log behavior of Ai^ 




Non-log behavior of Ai^ 


I 

III 
III' 
I' 


r < n < r 2 

i"i < r 2 < r 
r 2 < r\ < r 
r < r 2 < r\ 


icos(a)^ + 0(4^) 
ir 2 + 0(r 1 r 2 log(^)) 
ir 2 + 0(r 1 r 2 log(a?)) 
icos(a)^ + 0(4^) 

8 v ' r\ri ^r^r^' 


^r 4 

32 

^r 4 

32 ' 


^cos(a)^ + 0(4^) 
+ r 2 (0(r 2 ) + 0{nr 2 ) + 0(r 2 )) 
+ r 2 (0(r 2 ) + 0{ ri r 2 ) + 0(r 2 )) 

_L cos ( a )^ + 0(4^) 

24 V / riT2 \ r * r ^) 



References 

[1] P. F. Kolb, J. Sollfrank and U. W. Heinz, "Anisotropic transverse flow and the quark-hadron phase 
transition," Phys. Rev. C 62, 054909 (2000) [arXiv:hep-ph/0006129]. 

[2] P. F. Kolb, P. Huovinen, U. W. Heinz and H. Heiselberg, "Elliptic flow at SPS and RHIC: From 
kinetic transport to hydrodynamics," Phys. Lett. B 500, 232 (2001) [arXiv:hep-ph/0012137]. 

[3] P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, "Radial and elliptic flow 
at RHIC: further predictions," Phys. Lett. B 503, 58 (2001) [arXiv:hep-ph/0101136]. 



- 49 - 



[4] P. F. Kolb, U. W. Heinz, P. Huovinen, K. J. Eskola and K. Tuominen, "Centrality dependence of 
multiplicity, transverse energy, and elliptic flow from hydrodynamics," Nucl. Phys. A 696, 197 
(2001) [arXiv:hep-ph/0103234]. 

[5] U. W. Heinz and P. F. Kolb, "Early thermalization at RHIC," Nucl. Phys. A 702, 269 (2002) 
[arXiv:hep-ph/0111075]. 

[6] D. Teaney and E. V. Shuryak, "An unusual space-time evolution for heavy ion collisions at high 
energies due to the QCD phase transition," Phys. Rev. Lett. 83, 4951 (1999) 
[arXiv:nucl-th /9904006] . 

[7] D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. Lett. 86, 4783 (2001) [arXiv:nucl-th/0011058]. 

[8] D. Teaney, J. Lauret and E. V. Shuryak, "A hydrodynamic description of heavy ion collisions at the 
SPS and RHIC," arXiv:nucl-th/0110037. 

[9] D. Teaney, "Effect of shear viscosity on spectra, elliptic flow, and Hanbury Brown-Twiss radii," 
Phys. Rev. C 68, 034913 (2003) [arXiv:nucl-th/0301099]. 

[10] G. Policastro, D. T. Son and A. O. Starinets, "The shear viscosity of strongly coupled N = 4 
supersymmetric Yang-Mills plasma," Phys. Rev. Lett. 87, 081601 (2001) [arXiv:hep-th/0104066]. 

[11] D. T. Son and A. O. Starinets, "Minkowski-space correlators in AdS/CFT correspondence: Recipe 
and applications," JHEP 0209, 042 (2002) [arXiv:hep-th/0205051]. 

[12] G. Policastro, D. T. Son and A. O. Starinets, "From AdS/CFT correspondence to hydrodynamics," 
JHEP 0209, 043 (2002) [arXiv:hep-th/0205052]. 

[13] P. Kovtun, D. T. Son and A. O. Starinets, "Holography and hydrodynamics: Diffusion on stretched 
horizons," JHEP 0310, 064 (2003) [arXiv:hep-th/0309213]. 

[14] P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting quantum field theories 
from black hole physics," Phys. Rev. Lett. 94, 111601 (2005) [arXiv:hep-th/0405231]. 

[15] Y. V. Kovchegov and A. Taliotis, "Early time dynamics in heavy ion collisions from AdS/CFT 
correspondence," Phys. Rev. C 76, 014905 (2007) [arXiv:0705.1234 [hep-ph]]. 

[16] G. D. Moore and D. Teaney, "How much do heavy quarks thermalize in a heavy ion collision?," 
Phys. Rev. C 71, 064904 (2005) [arXiv:hep-ph/0412346]. 

[17] J. Casalderrey-Solana, E. V. Shuryak and D. Teaney, "Conical flow induced by quenched QCD jets," 
J. Phys. Conf. Ser. 27, 22 (2005) [Nucl. Phys. A 774, 577 (2006)] [arXiv:hep-ph/0411315]. 

[18] J. Casalderrey-Solana, E. V. Shuryak and D. Teaney, "Hydrodynamic flow from fast particles," 
arXiv:hep-ph/0602183. 

[19] E. V. Shuryak, "Strongly coupled quark-gluon plasma: The status report," arXiv:hep-ph/0608177. 



- 50 - 



R. Baier, A. H. Mueller, D. Scruff and D. T. Son, "'Bottom-up' thermalization in heavy ion 
collisions," Phys. Lett. B 502, 51 (2001) [arXiv:hep-ph/0009237]. 

P. B. Arnold, J. Lenaghan and G. D. Moore, "QCD plasma instabilities and bottom-up 
thermalization," JHEP 0308, 002 (2003) [arXiv:hep-ph/0307325]. 

S. Mrowczynski, "STREAM INSTABILITIES OF THE QUARK - GLUON PLASMA," Phys. Lett. 
B 214, 587 (1988) [Erratum-ibid. B 656, 273 (2007)]. 

S. Mrowczynski, Plasma instability at the initial stage of ultrarelativistic heavy ion collisions, Phys. 
Lett. B314 (1993) 118-121. 

Y. V. Kovchegov, Thoughts on non-perturbative thermalization and jet quenching in heavy ion 



collisions, Nucl. Phys. A764 (2006) 476-497, jhep-ph/0507134 l 



Y. V. Kovchegov, Isotropization and thermalization in heavy ion collisions, Nucl. Phys. A774 (2006) 
869-872, [ |hep-ph/0510232|1 . 



Y. V. Kovchegov, Can thermalization in heavy ion collisions be described by QCD diagrams?, Nucl. 



Phys. A762 (2005) 298-325, [hep-ph/0503038l. 



J. P. Blaizot and A. H. Mueller, The Early Stage of Ultrarelativistic Heavy Ion Collisions, Nucl. 
Phys. B289 (1987) 847. 

L. D. McLerran and R. Venugopalan, Computing quark and gluon distribution functions for very 



large nuclei, Phys. Rev. D49 (1994) 2233-2241, [)hep-ph/9309289 



L. D. McLerran and R. Venugopalan, Gluon distribution functions for very large nuclei at small 



transverse momentum, Phys. Rev. D49 (1994) 3352-3355, jhep-ph/9311205 



L. D. McLerran and R. Venugopalan, Green's functions in the color field of a large nucleus, Phys. 
Rev. D50 (1994) 2225-2233, [|hep-ph/9402335|1 . 



Y. V. Kovchegov, Non-abelian Weizsaecker- Williams field and a two- dimensional effective color 



charge density for a very large nucleus, Phys. Rev. D54 (1996) 5463-5469, [hep-ph/9605446|. 

Y. V. Kovchegov, Quantum structure of the non-abelian Weizsaecker-Williams field for a very large 



nucleus, Phys. Rev. D55 (1997) 5445-5455, |hep-ph/9701229 | 



A. Kovner, L. D. McLerran, and H. Weigert, Gluon production from nonAbelian 
Weizsacker-Williams fields in nucleus-nucleus collisions, Phys. Rev. D52 (1995) 6231-6237, 
Hhep-pn/9502289|1 . 

A. Krasnitz and R. Venugopalan, Non-perturbative computation of gluon mini-jet production in 
nuclear collisions at very high energies, Nucl. Phys. B557 (1999) 237, [|hep-ph/9809433|1 . 



A. Krasnitz and R. Venugopalan, The initial energy density of gluons produced in very high energy 



nuclear collisions, Phys. Rev. Lett. 84 (2000) 4309-4312, jhep-pn/9909203 1. 



- 51 - 



[36] 
[37] 
[38] 
[39] 
[40] 
[41] 

[42] 

[43] 
[44] 
[45] 
[46] 
[47] 
[48] 
[49] 
[50] 



A. Krasnitz, Y. Nara, and R. Venugopalan, Classical gluodynamics of high energy nuclear collisions: 



An erratum and an update, Nucl. Phys. A727 (2003) 427-436, [tiep-ph/0305112]. 



Y. V. Kovchegov, Classical initial conditions for ultrarelativistic heavy ion collisions, Nucl. Phys. 
A692 (2001) 557-582, [ frep-ph/00li~25^ . 



D. Kharzeev and M. Nardi, Hadron production in nuclear collisions at RHIC and high density QCD, 



Phys. Lett. B507 (2001) 121-128, [Emcl-th/00 12025 



D. Kharzeev, E. Levin, and M. Nardi, The onset of classical QCD dynamics in relativistic heavy ion 



collisions, Phys. Rev. C71 (2005) 054903, [hep-ph/0111315 



D. Kharzeev, E. Levin, and L. McLerran, Parton saturation and N(part) scaling of semi-hard 



processes in QCD, Phys. Lett. B561 (2003) 93-101, [hep-ph/0210332 



D. Kharzeev, Y. V. Kovchegov, and K. Tuchin, Nuclear modification factor in d + Au collisions: 
Onset of suppression in the color glass condensate, Phys. Lett. B599 (2004) 23-31, 
[)hep-pn/0405045|1 . 



J. L. Albacete, N. Armesto, A. Kovner, C. A. Salgado, and U. A. Wiedemann, Energy dependence of 
the Cronin effect from non-linear QCD evolution, Phys. Rev. Lett. 92 (2004) 082001, 
| |hep-ph/0307179| |. 

J. L. Albacete, Particle multiplicities in Lead-Lead collisions at the LHC from non-linear evolution 



with running coupling, Phys. Rev. Lett. 99 (2007) 262301, [|arXiv : 0707 . 2545 



E. Iancu and R. Venugopalan, The color glass condensate and high energy scattering in QCD, 
|hep-ph/0303204 . 



H. Weigert, Evolution at small x^: The Color Glass Condensate, Prog. Part. Nucl. Phys. 55 (2005) 



461-565, [tiep-ph/0501087]. 



J. Jalilian-Marian and Y. V. Kovchegov, Saturation physics and deuteron gold collisions at RHIC, 
Prog. Part. Nucl. Phys. 56 (2006) 104-231, Hhep-ph/0505052|1 . 



A. Krasnitz, Y. Nara, and R. Venugopalan, Gluon production in the color glass condensate model of 
collisions of ultrarelativistic finite nuclei, Nucl. Phys. A717 (2003) 268-290, |hep-ph/0209269 |. 



J. M. Maldacena, The large n limit of super conformal field theories and supergravity, Adv. Theor. 
Math. Phys. 2 (1998) 231-252, [|hep-th/97iT200|. 



E. Witten, Anti-de sitter space and holography, Adv. Theor. Math. Phys. 2 (1998) 253-291 
Hhep-th/9802150|1 . 



Y. Burnier, M. Laine and M. Vepsalainen, "Dimensionally regularized Polyakov loop correlators in 
hot QCD," JHEP 1001, 054 (2010) [arXiv:0911.3480 [hep-ph]]. 



- 52 - 



[51] J. Noronha and A. Dumitru, "The Heavy Quark Potential as a Function of Shear Viscosity at 
Strong Coupling," Phys. Rev. D 80, 014007 (2009) [arXiv:0903.2804 [hep-ph]]. 

[52] J. M. Maldacena, "Wilson loops in large N field theories," Phys. Rev. Lett. 80, 4859 (1998) 
[arXiv:hep-th/9803002]. 

[53] D. Bak, A. Karch and L. G. Yaffe, "Debye screening in strongly coupled N=4 supersymmetric 
Yang-Mills plasma," JHEP 0708, 049 (2007) [arXiv:0705.0994 [hep-th]]. 

[54] S. J. Rey, S. Theisen and J. T. Yee, "Wilson-Polyakov loop at finite temperature in large N gauge 
theory and anti-de Sitter supergravity," Nucl. Phys. B 527, 171 (1998) [arXiv:hep-th/9803135]. 

[55] J. L. Albacete, "Heavy Quark Potential at Finite Temperature in AdS/CFT," Nucl. Phys. A 830, 
311C (2009) [arXiv:0908.2541 [hep-ph]]. 

[56] J. L. Albacete, Y. V. Kovchegov and A. Taliotis, "Heavy Quark Potential at Finite Temperature in 
AdS/CFT Revisited," Phys. Rev. D 78, 115007 (2008) [arXiv:0807.4747 [hep-th]]. 

[57] A. Brandhuber, N. Itzhaki, J. Sonnenschein and S. Yankielowicz, "Wilson loops in the large N limit 
at finite temperature," Phys. Lett. B 434, 36 (1998) [arXiv:hep-th/9803137]. 

[58] C. Marquet, C. Roiesnel and S. Wallon, "Virtual Compton Scattering off a Spinless Target in 
AdS/QCD," arXiv: 1002.0566 [hep-ph]. 

[59] M. A. Betemps, V. P. Goncalves and J. T. de Santana Amaral, Phys. Rev. D 81, 094012 (2010) 
[arXiv: 1001. 3548 [hep-ph]]. 

[60] L. Cornalba and M. S. Costa, "Saturation in Deep Inelastic Scattering from AdS/CFT," Phys. Rev. 
D 78, 096010 (2008) [arXiv:0804.1562 [hep-ph]]. 

[61] L. Cornalba, M. S. Costa and J. Penedones, "Eikonal Methods in AdS/CFT: BFKL Pomeron at 
Weak Coupling," JHEP 0806, 048 (2008) [arXiv: 080 1.3002 [hep-th]]. 

[62] L. Cornalba, M. S. Costa and J. Penedones, "Deep Inelastic Scattering in Conformal QCD," 
arXiv:091 1.0043 [hep-th]. 

[63] J. Polchinski and M. J. Strassler, "Hard scattering and gauge / string duality," Phys. Rev. Lett. 88, 
031601 (2002) [arXiv:hep-th/0109174]. 

[64] J. Polchinski and M. J. Strassler, "Deep inelastic scattering and gauge/string duality," Phys. Rev. 
Lett. 88, 031601 (2002); JHEP 0305, 012 (2003). [arXiv:hep-th/0209211]. 

[65] J. H. Gao and B. W. Xiao, "Non-forward Compton scattering in AdS/CFT," Phys. Rev. D 81, 
035008 (2010) [arXiv:0912.4333 [hep-ph]]. 

[66] E. Iancu and A. H. Mueller, "Light-like mesons and deep inelastic scattering in finite-temperature 
AdS/CFT with flavor," JHEP 1002, 023 (2010) [arXiv:0912.2238 [hep-th]]. 



- 53 - 



[67] F. Dominguez, "Particle production in DIS off a in AdS," arXiv:0912.1641 [hep-th]. 

[68] Y. V. Kovchegov, "R-Current DIS on a Shock Wave: Beyond the Eikonal Approximation," 
arXiv:1005.0374 [hep-ph]. 

[69] J. L. Albacete, Y. V. Kovchegov and A. Taliotis, "DIS on a Large Nucleus in AdS/CFT," JHEP 
0807, 074 (2008) [arXiv:0806.1484 [hep-th]]. 

[70] J. L. Albacete, Y. V. Kovchegov and A. Taliotis, "DIS in Ads," AIP Conf. Proc. 1105, 356 (2009) 
[arXiv:0811.0818 [hep-th]]. 

[71] A. H. Mueller, A. I. Shoshi and B. W. Xiao, "Deep inelastic and dipole scattering on finite length 
hot N = 4 SYM matter," Nucl. Phys. A 822, 20 (2009) [arXiv:0812.2897 [hep-th]]. 

[72] C. Marquet and T. Renk, "Jet quenching in the strongly-interacting quark-gluon plasma," Phys. 
Lett. B 685, 270 (2010) [arXiv:0908.0880 [hep-ph]]. 

[73] C. A. Ballon Bayona, H. Boschi-Filho and N. R. F. Braga, "Deep inelastic scattering from gauge 
string duality in D3-D7 brane model," JHEP 0809, 114 (2008) [arXiv:0807.1917 [hep-th]]. 

[74] A. Taliotis, "DIS from the AdS/CFT correspondence," Nucl. Phys. A 830, 299C (2009) 
[arXiv:0907.4204 [hep-th]]. 

[75] M. Giordano and R. Peschanski, "High Energy Bounds on Soft N=4 SYM Amplitudes from 
AdS/CFT," arXiv:1003.2309 [hep-ph]. 

[76] Y. Hatta, "Relating e + e~ annihilation to high energy scattering at weak and strong coupling," 
JHEP 0811, 057 (2008) [arXiv:0810.0889 [hep-ph]]. 

[77] E. Avsar, E. Iancu, L. McLerran and D. N. Triantafyllopoulos, "shock waves and deep inelastic 
scattering within the gauge/gravity duality," JHEP 0911, 105 (2009) [arXiv:0907.4604 [hep-th]]. 

[78] L. Cornalba, M. S. Costa and J. Penedones, "AdS black disk model for small-x DIS," 
arXiv: 1001. 1157 [hep-ph]. 

[79] R. C. Brower, M. Djuric, I. Sarcevic and C. I. Tan, "String-Gauge Dual Description of Deep 
Inelastic Scattering at Small-x," arXiv:1007.2259 [hep-ph]. 

[80] Y. V. Kovchegov, Z. Lu and A. H. Rezaeian, "Comparing AdS/CFT Calculations to HERA F 2 
Data," Phys. Rev. D 80, 074023 (2009) [arXiv:0906.4197 [hep-ph]]. 

[81] W. A. Horowitz and Y. V. Kovchegov, "Shock Treatment: Heavy Quark Drag in a Novel AdS 
Geometry," Phys. Lett. B 680, 56 (2009) [arXiv:0904.2536 [hep-th]]. 

[82] W. A. Horowitz and M. Gyulassy, "Heavy quark jet tomography of Pb + Pb at LHC: AdS/CFT 
drag or pQCD energy loss?," Phys. Lett. B 666, 320 (2008) [arXiv: 0706. 2336 [nucl-th]]. 



-54- 



W. A. Horowitz, "Shock Treatment: Heavy Quark Energy Loss in a Novel AdS/CFT Geometry," 
Nucl. Phys. A 830, 773C (2009) [arXiv:0907.4845 [nucl-th]]. 

S. S. Gubser, D. R. Gulotta, S. S. Pufu and F. D. Rocha, "Gluon energy loss in the gauge-string 
duality," JHEP 0810, 052 (2008) [arXiv:0803.1470 [hep-th]]. 

S. D. Avramis, K. Sfetsos and D. Zoakos, "On the velocity and chemical-potential dependence of the 
heavy-quark interaction in N = 4 SYM plasmas," Phys. Rev. D 75, 025009 (2007) 
[arXiv:hep-th/0609079]. 

C. Athanasiou, P. M. Chesler, H. Liu, D. Nickel and K. Rajagopal, "Synchrotron radiation in 
strongly coupled conformal field theories," arXiv:1001.3880 [hep-th]. 

G. Beuf, C. Marquet and B. W. Xiao, "Heavy-quark energy loss and thermalization in a strongly 
coupled SYM plasma," Phys. Rev. D 80, 085001 (2009) [arXiv:0812.1051 [hep-ph]]. 

K. Kajantie, T. Tahkokallio, and J.-T. Yee, Thermodynamics of ads/qcd, JHEP 01 (2007) 019, 
[)hep-pri/0609254| l. 

K. Kajantie and T. Tahkokallio, Spherically expanding matter in ads/cft, Phys. Rev. D75 (2007) 
066003, flriep-th/0612226 ]. 

K. Kajantie, J. Louko, and T. Tahkokallio, The gravity dual of 1+1 dimensional bjorken expansion, 
|0705.179i . 



K. Kajantie, J. Louko, and T. Tahkokallio, Gravity dual of conformal matter collisions in 1+1 
dimensions, Phys. Rev. D77 (2008) 066001, [|arXiv:0801 .0198| 1. 

J. L. Albacete, Y. V. Kovchegov and A. Taliotis, "Asymmetric Collision of Two Shock Waves in 
AdS 5 ," JHEP 0905, 060 (2009) [arXiv:0902.3046 [hep-th]]. 

J. L. Albacete, Y. V. Kovchegov and A. Taliotis, "Modeling Heavy Ion Collisions in AdS/CFT," 
JHEP 0807, 100 (2008) [arXiv:0805.2927 [hep-th]]. 

D. Grumiller and P. Romatschke, "On the collision of two shock waves in AdS5," JHEP 0808, 027 
(2008) [arXiv:0803.3226 [hep-th]]. 

R. A. Janik and R. Peschanski, Asymptotic perfect fluid dynamics as a consequence of AdS/CFT, 
Phys. Rev. D73 (2006) 045013, [|hep-th/0512162|1 . 

R. A. Janik, "The dynamics of quark-gluon plasma and AdS/CFT," arXiv: 1003.3291 [hep-th]. 

R. A. Janik and R. Peschanski, Gauge / gravity duality and thermalization of a boost- invariant 
perfect fluid, Phys. Rev. D74 (2006) 046007, piep-th/0606149 |. 

[98] M. P. Heller, R. A. Janik and R. Peschanski, "Hydrodynamic Flow of the Quark-Gluon Plasma and 
Gauge/Gravity Correspondence," Acta Phys. Polon. B 39, 3183 (2008) [arXiv:0811.3113 [hep-th]]. 



- 55 - 



[99] G. Beuf, M. P. Heller, R. A. Janik, and R. Peschanski, Boost-invariant early time dynamics from 
[100] 
[101] 



[102 

[103 

[104 

[105 

[106 

[107 

[108 
[109 

[110 

[111 

[112 

[113 
[114 



AdS/CFT, JEEP 10 (2009) 043, [|arXiv: 0906. 4423 ]. 

R. A. Janik, Viscous plasma evolution from gravity using AdS/CFT, Phys. Rev. Lett. 98 (2007) 
022302, jhep-th/0610144 1. 



D. Bak and R. A. Janik, From static to evolving geometries: R-charged hydrodynamics from 
supergravity, Phys. Lett. B645 (2007) 303-308, [ |iep-tri/0611304| ]. 

S. Nakamura and S.-J. Sin, A holographic dual of hydrodynamics, JEEP 09 (2006) 020, 
Hhep-th/060712311 . 



E. Shuryak, S.-J. Sin, and I. Zahed, A gravity dual of REJC collisions, J. Korean Phys. Soc. 50 
(2007) 384-397, | fcep-th/0511199| |. 

S. Lin and E. Shuryak, Toward the AdS/CFT gravity dual for high energy heavy ion collisions, 
|hep-ph/0610168j . 



M. Lublinsky and E. Shuryak, Eow much entropy is produced in strongly coupled quark- gluon 
plasma (sQGP) by dissipative effects?, p704 . 1647 . 

I. Y. Aref'eva, "Catalysis of Black Holes/ Wormholes Formation in High Energy Collisions," 
arXiv:0912.5481 [hep-th]. 

O. Aharony, S. Minwalla, and T. Wiseman, Plasma-balls in large N gauge theories and localized 
black holes, Class. Quant. Grav. 23 (2006) 2171-2210, jhep-th/0507219|] . 

S. S. Gubser, arXiv:1006.0006 [hep-th]. 

S. Bhattacharyya, R. Loganayagam, S. Minwalla, S. Nampuri, S. P. Trivedi and S. R. Wadia, 
"Forced Fluid Dynamics from Gravity," JHEP 0902, 018 (2009) [arXiv:0806.0006 [hep-th]]. 

T. Springer, "Second order hydrodynamics for a special class of gravity duals," Phys. Rev. D 79, 
086003 (2009) [arXiv: 0902. 2566 [hep-th]]. 

T. Springer, "Sound Mode Hydrodynamics from Bulk Scalar Fields," Phys. Rev. D 79, 046003 
(2009) [arXiv:0810.4354 [hep-th]]. 

J. I. Kapusta and T. Springer, "Shear Transport Coefficients from Gauge/Gravity 
Correspondence," Phys. Rev. D 78, 066017 (2008) [arXiv:0806.4175 [hep-th]]. 

H. Nastase, The REJC fireball as a dual black hole, [hep-th/ 0501068] . 

S. S. Gubser, S. S. Pufu, and A. Yarom, Entropy production in collisions of gravitational shock 
waves and of heavy ions, Phys. Rev. D78 (2008) 066014, [ |arXiv : 0805 . 155l| . 



[115] S. S. Gubser, S. S. Pufu, and A. Yarom, Off-center collisions in AdS§ with applications to 
multiplicity estimates in heavy-ion collisions, JEEP 11 (2009) 050, | arXiv:0902 .4062]| . 



- 56 - 



[116 

[ii7; 
[us; 

[119 

[120 
[121 
[122 
[123 
[124 

[125; 

[126 
[127 

[128; 

[129 
[130 



S. Lin and E. Shuryak, Grazing Collisions of Gravitational Shock Waves and Entropy Production in 



Heavy Ion Collision, Phys. Rev. D79 (2009) 124015, [arXiv: 0902 . 1508]. 



Y. V. Kovchegov and S. Lin, "Toward Thermalization in Heavy Ion Collisions at Strong Coupling," 
JHEP 1003, 057 (2010) [arXiv:0911.4707 [hep-th]]. 

L. Alvarez-Gaume, C. Gomez, A. Sabio Vera, A. Tavanfar, and M. A. Vazquez-Mozo, Critical 
formation of trapped surfaces in the collision of gravitational shock waves, JHEP 02 (2009) 009, 



arXiv:0811.3969 



I. Y. Aref'eva, A. A. Bagrov and E. A. Guseva, "Critical Formation of Trapped Surfaces in the 
Collision of Non-expanding Gravitational Shock Waves in de Sitter Space-Time," JHEP 0912, 009 
(2009) [arXiv:0905.1087 [hep-th]]. 

I. Y. Aref'eva, A. A. Bagrov and L. V. Joukovskaya, "Critical Trapped Surfaces Formation in the 
Collision of Ultrarelativistic Charges in (A)dS," JHEP 1003, 002 (2010) [arXiv:0909.1294 [hep-th]]. 

P. M. Chesler and L. G. Yaffe, Horizon formation and far- from- equilibrium isotropization in 



super symmetric Yang- Mills plasma, Phys. Rev. Lett. 102 (2009) 211601, |arXiv:0812.2053l. 



P. M. Chesler and L. G. Yaffe, Boost invariant flow, black hole formation, and far-from- equilibrium 
dynamics in N = 4 super symmetric Yang-Mills theory, |arXiv : 0906 . 4426 . 



A. J. Amsel, D. Marolf, and A. Virmani, "Collisions with Black Holes and Deconfined Plasmas," 



JHEP 04 (2008) 025, |0712.2221 



H. Nastase, AdS-CFT and the RHIC fireball, Prog. Theor. Phys. Suppl. 17 A (2008) 274-285, 
flarXiv: 0805 .3579 ]. 



S. B. Giddings and R. A. Porto, "The gravitational S-matrix," Phys. Rev. D 81, 025002 (2010) 
[arXiv:0908.0004 [hep-th]]. 

S. B. Giddings, M. Schmidt-Sommerfeld and J. R. Andersen, "High energy scattering in gravity 
and super gravity," arXiv: 1005.5408 [hep-th]. 

D. M. Eardley and S. B. Giddings, "Classical black hole production in high-energy collisions," 
Phys. Rev. D 66, 044011 (2002) [arXiv:gr-qc/0201034]. 

S. B. Giddings and V. S. Rychkov, "Black holes from colliding wavepackets," Phys. Rev. D70 
(2004) 104026, [hep-th/040"9T3"H . 



M. Panero, "Thermodynamics of the QCD plasma and the large-N limit," Phys. Rev. Lett. 103, 
232001 (2009) [arXiv:0907.3719 [hep-lat]]. 

U. Giirsoy and E. Kiritsis, Exploring improved holographic theories for QCD: Part I, JHEP 0802, 
032 (2008) [arXiv:0707.1324 [hep-th]]. 



-57- 



[131 
[132 
[133 
[134 
[135 
[136 

[137; 

[138 

[139 
[140 
[141 
[142 
[143 
[144 
[145 



U. Giirsoy, E. Kiritsis and F. Nitti, Exploring improved holographic theories for QCD: Part II, 
JHEP 0802, 019 (2008) [arXiv:0707.1349 [hep-th]]. 

U. Giirsoy, E. Kiritsis, L. Mazzanti and F. Nitti, Deconfinement and gluon plasma dynamics in 
improved holographic QCD, Phys. Rev. Lett. 101, 181601 (2008) [arXiv:0804.0899 [hep-th]]. 

E. Kiritsis, Dissecting the string theory dual of QCD, Fortsch. Phys. 57, 396 (2009) 
[arXiv:0901.1772 [hep-th]]. 

J. D. Edelstein, J. P. Shock and D. Zoakos, "The AdS/CFT Correspondence and Non-perturbative 
QCD," AIP Conf. Proc. 1116, 265 (2009) [arXiv: 090 1.2534 [hep-ph]]. 

G. 't Hooft, "A PLANAR DIAGRAM THEORY FOR STRONG INTERACTIONS," Nucl. Phys. 
B 72, 461 (1974). 

S. de Haro, S. N. Solodukhin, and K. Skenderis, Holographic reconstruction of space-time and 
renormalization in the AdS/CFT correspondence, Commun. Math. Phys. 217 (2001) 595-622, 
[ piep-th/0002230| ]. 



P. D. D'Eath and P. N. Payne, "Gravitational radiation in high speed black hole collisions. 1. 
Perturbation treatment of the axisymmetric speed of light collision," Phys. Rev. D46 (1992) 
658-674. 

P. D. D'Eath and P. N. Payne, "Gravitational radiation in high speed black hole collisions. 2. 
Reduction to two independent variables and calculation of the second order news function," Phys. 
Rev. D46 (1992) 675-693. 

P. D. D'Eath and P. N. Payne, "Gravitational radiation in high speed black hole collisions. 3. 
Results and conclusions," Phys. Rev. D46 (1992) 694-701. 

P. C. Aichelburg and R. U. Sexl, "On the Gravitational field of a massless particle," Gen. Rel. 
Grav. 2, 303 (1971). 

T. Dray and G. 't Hooft, "The Gravitational Shock Wave of a Massless Particle," Nucl. Phys. 
B253 (1985) 173. 

M. Hotta and M. Tanaka, Shock wave geometry with nonvanishing cosmological constant, Class. 
Quant. Grav. 10 (1993) 307314. 

S. W. Hawking and R. Penrose, The Singularities of gravitational collapse and cosmology, Proc. 
Roy. Soc. Lond. A314 (1970) 529-548. 

S. Khlebnikov, M. Kruczenski and G. Michalogiorgakis, "Shock waves in strongly coupled 
plasmas," arXiv: 1004.3803 [hep-th]. 

K. Sfetsos, "On gravitational shock waves in curved space-times," Nucl. Phys. B436 (1995) 



721-746, |hep-th/9408169 . 



- 58 - 



[146 
[147 
[148 
[149 
[150 
[151 
[152 
[153 
[154 
[155 
[156 

[157; 
[158; 

[159 
[160 
[161 



G. Beuf, "Gravity dual of N=4 SYM theory with fast moving sources," Phys. Lett. B 686, 55 
(2010) [arXiv:0903.1047 [hep-th]]. 

T. Lappi, Production of gluons in the classical field model for heavy ion collisions, Phys. Rev. C67 
(2003) 054903, [[hep-ph/0303076 |. 



K. Fukushima, Initial fields and instability in the classical model of the heavy-ion collision, 
|0704.3625 . 



R. J. Fries, J. I. Kapusta, and Y. Li, Near-fields and initial energy density in the color glass 



condensate model, nucl-th/0604054 



T. Lappi and R. Venugopalan, "Universality of the saturation scale and the initial eccentricity in 
heavy ion collisions," Phys. Rev. C 74, 054905 (2006) [arXiv:nucl-th/0609021]. 

W. Jas and S. Mrowczynski, "Evolution of Anisotropy of Parton System from Relativistic 
Heavy-Ion Collisions," Phys. Rev. C 76, 044905 (2007) [arXiv:0706.2273 [nucl-th]]. 

P. F. Kolb and U. W. Heinz, "Hydrodynamic description of ultrarelativistic heavy-ion collisions," 
arXiv:nucl-th/0305084. 

P. F. Kolb and U. W. Heinz, "Emission angle dependent HBT at RHIC and beyond," Nucl. Phys. 
A 715, 653 (2003) [arXiv:nucl-th/0208047]. 

J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, "QCD and a holographic model of hadrons," 
Phys. Rev. Lett. 95, 261602 (2005) [arXiv:hep-ph/0501128]. 

A. Karch, E. Katz, D. T. Son and M. A. Stephanov, "Linear confinement and AdS/QCD," Phys. 
Rev. D 74, 015005 (2006) [arXiv:hep-ph/0602229]. 

J. Hirn, N. Rius and V. Sanz, "Geometric approach to condensates in holographic QCD," Phys. 
Rev. D 73, 085005 (2006) [arXiv:hep-ph/0512240]. 

J. Erlich, G. D. Kribs and I. Low, "Emerging holography," Phys. Rev. D 73, 096001 (2006). 
[arXiv:hep-th/0602110]. 

J. Hirn and V. Sanz, "Interpolating between low and high energy QCD via a 5D Yang-Mills 
model," JHEP 0512, 030 (2005). [arXiv:hep-ph/0507049]. 

L. Da Rold and A. Pomarol, "Chiral symmetry breaking from five dimensional spaces," Nucl. Phys. 
B 721, 79 (2005); [arXiv:hep-ph/0501218]. 

H. Boschi-Filho and N. R. F. Braga, "Gauge / string duality and scalar glueball mass ratios," 
JHEP 0305, 009 (2003); [arXiv:hep-th/0212207]. 

S. J. Brodsky and G. F. de Teramond, "Light-front hadron dynamics and AdS/CFT 
correspondence," Phys. Lett. B 582, 211 (2004); [apXiv:hep-th/0310227]. 



- 59 - 



[162 
[163 
[164 
[165 
[166 
[167; 

[168 
[169 
[170 
[171 
[172; 

[173; 

[174 
[175 



S. Hong, S. Yoon and M. J. Strassler, "On the couplings of vector mesons in AdS/QCD," JHEP 
0604, 003 (2006). [arXiv:hep-th/0409118]. 

T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," Prog. Theor. Phys. 
113, 843 (2005); [arXiv:hep-th/0412141]. 

T. Hambye, B. Hassanain, J. March-Russell and M. Schvellinger, "On the Delta(I) = 1/2 rule in 
holographic QCD," Phys. Rev. D 74, 026003 (2006). 

K. Ghoroku, N. Maru, M. Tachibana and M. Yahiro, "Holographic model for hadrons in deformed 
AdS(5) background," Phys. Lett. B 633, 602 (2006). [arXiv:hep-ph/0510334]. 

N. Evans, A. Tedder and T. Waterson, "Improving the infra-red of holographic descriptions of 
QCD," JHEP 0701, 058 (2007). [arXiv:hep-ph/0603249]. 

S. K. Domokos and J. A. Harvey, "Baryon number-induced Chern-Simons couplings of vector and 
axial- vector mesons in holographic QCD," Phys. Rev. Lett. 99, 141602 (2007). [arXiv:0704.1604 
[hep-ph]]. 

H. J. Kwee and R. F. Lebed, "Pion Form Factors in Holographic QCD," JHEP 0801, 027 (2008) 
[arXiv:0708.4054 [hep-ph]]. 

H. R. Grigoryan and A. V. Radyushkin, "Form Factors and Wave Functions of Vector Mesons in 
Holographic QCD," Phys. Lett. B 650, 421 (2007) [arXiv:hep-ph/0703069]. 

H. R. Grigoryan and A. V. Radyushkin, "Pion Form Factor in Chiral Limit of Hard- Wall 
AdS/QCD Model," Phys. Rev. D 76, 115007 (2007) [arXiv:0709.0500 [hep-ph]]. 

I. S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, Fifth Edition , 
(Academic Press, San Diego, 1994) , Eqs. (6.578.11) and (8.751.3). 

A. Papapetrou, " Lectures on general Relativity," (D. Reidel publishing Company, Holland 1974). 
U. H. Danielsson, E. Keski-Vakkuri, and M. Kruczenski, Vacua, Propagators, and Holographic 



Probes in AdS/CFT, JHEP 01 (1999) 002, 1hep-th/9812007|] . 



Y. V. Kovchegov and D. H. Rischke, "Classical gluon radiation in ultrarelativistic nucleus nucleus 
collisions," Phys. Rev. C 56, 1084 (1997) [arXiv:hep-ph/9704201]. 

A. Taliotis, "Evolving Geometries in General Relativity," In Partial Fulfillment of the Requirements 
for the Degree MASTER of Science in the Department of Mathematics of The Ohio State 
University, arXiv: 1007. 1452 [hep-th]. 



- 60 - 



